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ABSTRACT 


The  broad  objective  of  this  research  is  to  develop  techniques  to  reduce 
non-specular  scattering  from  objects  of  interest  to  the  Air  Force.  The  major 
tools  used  thus  far  in  the  investigations  have  been  two-dimensional  computer 
programs  which  solve  the  integral  equations  for  a  surface  impedance  boundary 
condition  and  for  a  thin  resistive  sheet  placed  near  a  conducting  obstacle.  The 
existing  programs  have  proven  valuable,  but  the  research  requirements  have 
now  outstripped  their  capabilities  and  more  inclusive  programs  are  required. 

The  main  Rirust  of  this  Interim  Report  is  upon  the  development  of  the  integral 
equations  necessary  for  the  construction  of  more  comprehensive  programs,  but 
for  the  sake  of  completeness,  the  equations  for  the  existing  programs  are  also 
derived  as  a  matter  of  course.  Of  considerable  importance  is  the  treatment  of 
"magnetic"  resistive  sheets  in  addition  to  conventional  electric  resistive  sheets 
so  that,  with  but  a  slight  expansion  of  the  concept  of  resistivity,  magneto-dielec¬ 
tric  layers  of  physical  materials  may  be  modeled.  In  addition  to  the  theoretical 
work,  experimental  results  obtained  bv  both  The  University  of  Michigan  Radiation 
Laboratory  and  the  Air  Force  Avionics  Laboratory  are  compared  with  computed 
predictions. 
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INTRODUCTION 

’  The  objective  of  the  research  described  in  this  Interim  Report  is  to  develop 

techniques  to  reduce  non-specuiar  scattering  from  typical  aerospace  vehicles  of 
interest  to  the  Air  Force,  Such  scattering  arises  from  both  smooth  and  edged 
structures  when  viewed  in  non-specuKr  directions  and  during  this  reporting 
period  an  ogival  cylinder  has  been  used  as  the  basic  shape  with  which  to  explore 
promising  techniques.  The  baric  tools  in  the  research  were  a  pair  of  computer 
programs  that  solve  the  appropriate  integral  equations  for  the  surface  fields  in 
i  two  dimensions. 

i 

i  The  programs  are  modified  versions  of  programs  furnished  by  AFAL,  the 

modifications  being  necessary  because  the  originals  lacked  the  flexibility  to  con¬ 
veniently  accommodate  the  shapes  of  interest  in  the  present  study.  One  of  the 
programs  (RAMD)  is  itself  an  outgrowth  of  an  cailier  modification  which  has 

t 

now  been  abandoned  (RAM IB)  and  is  based  upon  the  integral  equations  in  which 
a  surface  impedance  boundary  condition  is  imposed.  RAMD  solves  the  equations 
for  either  F  •  or  H-polarization  (incident  electric  or  magnetic  vector  parallel  to 
the  cylinder  axis,  respectively)  and  a  brief  description  of  it,  along  with  a  source 
listing,  is  given  in  ApDendix  B.  The  other  program  (REST)  sol'.es  the  integral 
equations  for  a  metallic  cylinder  which  may  have  resistioe  sheets  placed  on  or 
near  the  surface,  but  is  applicable  only  for  E-polarization.  A  source  listing  of 

! 

REST  may  be  found  in  Appendix  B  of  Knott,  Liepa  and  Senior  (1973). 

The  impetus  to  discard  the  now  defunct  RAM1B  developed  late  in  the 
j  predecessor  Contract  (F33615-72-C-1439)  when  it  became  apparent  that  the 

I  program  would  never  satisfactorily  produce  the  correct  solution  for  edged 

j  bodies  for  E-polarization.  At  the  time  the  Final  Report  for  that  Contract  was 

i 

I 

j  being  prepared,  we  were  still  trying  to  obtain  acceptable  values  by  sampling 

■  the  surface  densely  in  the  vicinity  of  the  edges.  Although  the  results  tended  to 

!  converge  on  what  was  believed  (and  has  been  since  demonstrated)  to  be  the  correct 

I 
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behavior,  the  sampling  rate  over  the  remainder  of  the  profile  had  to  be  reduced. 

A  summary  of  these  efforts  and  the  consequences  of  coarse  sampling  is  given  in 
Chapter  V;  not  unexpectedly,  accuracy  is  reduced  and  ultimately  the  gain  in 
accuracy  by  dense  packing  at  the  edges  will  be  offset  by  the  loss  due  to  sparse 
sampling  over  the  rest  of  the  body. 

The  shortcomings  of  RAM1B  were  serious  enough  to  demand  a  thorough 
examination  of  the  integral  equations  for  the  impedance  boundary  eordition,  and 
in  Chapter  II  the  equations  are  derived  using  both  the  scalar  and  vector  approaches. 
Because  the  impedance  boundary  condition  rigidly  links  the  electee  surface  currents 
to  the  magnetic  ones,  a  pair  of  integral  equations  is  produced  for  each  polarization 
and,  from  a  mathematical  viewpoint,  one  is  as  good  as  the  other.  From  a 
numerical  standpoint,  however,  one  of  the  pair  is  markedly  superior  and  al¬ 
though  RAM1B  used  the  "good''  equation  for  one  polarization,  it  used  the  "bad" 
one  for  the  other.  This  discovery  led  to  the  construction  of  program  RAMD, 
which  is  sufficiently  different  from  its  predecessor  to  be  regarded  as  a  new 
program,  and  not  merely  a  modification. 

Program  REST  has  been  the  basic  tool  in  a  study  still  in  progress  exploring 
the  use  of  resistive  sheets  for  non-specular  scattering  reduction,  the  results  of 
which  will  be  documented  elsewhere.  The  program  is  designed  only  for  E-polari¬ 
zation,  which  turns  out  to  be  the  case  when  a  simple  integral  equation  is  sufficient. 
The  H-polarization  integral  equation,  on  the  other  hand,  contains  a  second  derivative 
singularity  (and  is  thus  a  "bad"  equation)  that  poses  annoying,  but  surmountable, 
numerical  problems,  and  it  was  probably  for  this  reason  that  the  original  version 
(RAMC)  was  vestricted  to  E-polarization.  Since  both  polarizations  must  be 
treatable  in  order  to  carry  out  the  Contract  objectives,  we  undertook  the  theo¬ 
retical  study  of  Chapter  IV  to  secure  the  necessary  equations..  Unlike  the 
impedance  boundary  condition  case,  theie  is  no  choice;  only  one  equation  is 
prod  ced  for  each  polarization  be  it  simple  or  not,  and  it  must  be  used  to  obtain 
the  solution.  We  should  mention  in  passing  that  the  study  in  Chapter  in,  which 
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considers  impedance  boundary  condition  sheets,  has  provided  a  key  result  for 
the  analysis  and  interpretation  of  the  data  generated  by  program  REST. 

As  pointed  out  in  Chapter  IV,  the  principle  of  duality  permits  us  to 
immediately  write  down  the  equations  for  a  "magnetic"  resistance  sheet  from 
the  electric  sheet  equations.  Since  the  magnetic  and  electric  currents  are 
independent  for  each  kind  of  sheet,  the  sheets  may  be  mathematically  super¬ 
posed  to  simulate  a  thin  layer  of  magneto-dielectric  material,  and  we  now  have 
the  foundation  for  a  computer  program  capable  of  modeling  physicaJy  realizable 
materials  in  contrast  to  the  more  nebulous  specification  of  surface  impedances. 
The  numerical  evaluation  of  second  order  singularities  is  an  inevitable  require¬ 
ment  in  a  generalized  program,  of  course,  but  the  procedures  discussed  in 
Appendix  A  show  that  the  task  can  be  accomplished  with  a  tolerably  small  error. 
And  although  the  thickness  of  actual  or  desirable  coatings  may  well  prove  to  be 
beyond  the  bounds  of  a  thin  layer  approximation,  a  thick  coating  can  be  mathe¬ 
matically  modeled  by  a  stack  of  thin  layers. 

The  need  frr  such  a  generalized  program  may  become  apparent  from  the 
compai'isons  between  measured  and  computed  results  given  in  Chapter  V.  The 
computations  were  performed  by  program  RAMD  and  although  they  agree  well 
with  experiment  for  metallic,  and  in  some  cases  coated,  bodies,  the  agreement 
for  other  coatings  is  not  as  good  The  disparity  in  these  cases  is  not  necessarily 
the  fault  of  the  program,  nor  indeed  of  the  experiments,  but  a  consequence  of  the 
difficulty  of  correctly  relating  the  impedance  near  an  edge  to  the  bulk  properties 
of  a  material  coating  and  the  local  edge  geometry.  This  specification  in  itself 
requires  the  solution  oi  a  boundary  value  problem,  a  task  we  hope  to  accomplish 
via  a  generalized  computer  program  based  upon  the  integral  equations  developed 
in  this  report. 
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INTEGRAL  EQUATION.'-  FOR  IMPEDANCE  SURFACES 

It  is  desirable  to  examine  the  form  of  integral  equations  appropriate  in 
a  two-dimensional  problem  involving  a  scatterer  with  an  impedance  boundary 
condition  imposed  at  its  surface.  The  objectives  of  this  study  are  not  only  to 
show  how  the  equations  are  affected  by  the  geometry  and  the  boundary  conditions, 
but  also  to  resolve  some  of  the  difficulties  which  have  previously  been  experienced. 

2. 1  The  Impedance  Boundary  Condition 

Consider  a  closed  cylindrical  surface  whose  generators  are  parallel  to 
the  z  axis  of  a  cylindrical  coordinate  system  and  whose  profile  in  a  plane  perpen¬ 
dicular  to  the  z  axis  is  the  closed  curve  C.  Let  n  be  a  unit  vector  normal  drawn 
outwards  to  C  and  let  s  be  a  unit  vector  in  the  tangential  direction,  such  that 
s,  n,  z  form  a  right-handed  system. 

At  the  surface  the  impedance  boundary  condition 

E-(n-E)n  =r]Z0nAH  (2.1) 

is  imposed,  where  rj  is  the  (normalized)  surface  impedance  and  may  be  a  function 
of  the  distance  s  along  C..  For  a  perfectly  conducting  surface,  n  =  0.  By  trivial 
manipulation,  eq.  (2, 1)  can  be  expressed  alternatively  as 

Y 

H  -  (n  •  H)  n  =  -  —  n  .  E  (2.2) 

-  —  n  A  — 

(Senior,  1962),  which  is  the  "dual"  of  (2.  1)  under  the  transformation  E/~>H_, 

II— >-E  ,  — >  Yq,  ri  < — >  l/r]..  Of  course,  this  duality  is  lost  if  at  any  p^int  of 

the  surface  t)  is  either  zero  or  infinite. 

The  body  is  illuminated  by  a  plane  wave  whoso  direction  of  propagation  is 
in  a  plane  perpendicular  to  the  z  axis,  and  it  is  now  sufficient  to  consider  separ¬ 
ately  the  two  principal  polarizations,  that  in  which  H*  =  zH*  (H-polaiization), 
and  E1  =  zE*  (E-polarization).  From  symmetry  it  follows  that  the  total  fields 
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are  similarly  oriented  and  this  enables  us  to  treat  each  problem  as  a  scalar  one 
for  the  z  component  of  the  total  magnetic  or  electric  field.  Such  a  simplification 
is  not  possible  if  the  plane  wave  is  at  skew  incidence. 

2.2  Scalar  Formulation 

Let  us  consider  first  the  case  of  H-polarization.  Then  H1  =  zH1, 

z 

implying  if  =  z  H  ,  and  consonant  with  the  usual  definition  of  the  surface  electric 
current  K,  viz. 

K=nAH,  (2,3) 

we  have 

K  =  sK  with  K  =  H  .  (2.4) 

—  s  s  z 

On  the  surface  the  boundary  condition  (2.  1)  requires 


E  =  rjZ  H 
s  0  z 


bZ0K 


s 


and  since,  from  Maxwell’s  equations. 


=  -ikY.E 
0  s 


we  have 


9H 

=  -ikr?K  on  C.  (2.  5) 

9n  s 

Application  of  the  scalar  Green’s  theorem  to  a  function  </  (g)  satisfying  the 
scalar  wave  equation  outside  and  on  the  surface  gives 


</(£) 


*l<e> 


fd1)(kr)-Hi1)(kr)^7^(g’) 


0 


*0 


ds* 


(2.6) 


J  C  J 

where  the  normal  derivatives  are  with  respect  to  the  primed  coordinates  of  the 
integration  point,  £  is  the  two-dimensional  position  vector  of  the  observation  point 


o 


ana  r  =  implying  r  =  j£-£'  | ..  For  simplicity  we  shall  henceforth  specify 

the  integration  point  by  its  circumferential  distance  s'  along  C. 

If  we  identify  </•  with  H  ,  eq.  (2.  6)  becomes 


r)H^(kr)+in(s')  H^’&r) 


K  (s’)d(ks*) 

s 


on  making  use  of  (2.  5),  and  if  we  now  allow  £  to  approach  a  point  s  on  C,  the 
following  integral  equation  results: 


H1  (s)  =  K  (s)  + 
z  s  £- 


r}(s*)H^(kr)-i(n'  •  r)H^(kr) 


K  (s')d(ks1). 
s 


The  second  term  in  the  integrand  has  a  non-integrable  singularity  when  £  is  on 
C,  but  if  this  is  treated  analytically  in  the  usual  manner,  the  limit  can  be  applied 
to  the  integrand  to  yield 


H^s)  =  ~  K  (s)  +  7 
z  2  s  4 


\  rj(s')K  (sOH^Ckrldlks') 
»  s  0 

Jc 


-7  X  K  (s'Hn*  •  r)H^1)(kr)d(ksl)  .  (2.9) 

4  \  s  1 

JC 

where  the  slash  across  the  integral  sign  indicates  a  Cauchy  principal  value.  This 
is  identical  to  the  equation  obtained  by  Andreasen  (1965)  for  the  special  case  of 
incidence  in  a  plane  perpendicular  to  the  z  axis  and  u as  quoted  (eq.  3.  l)  in 
the  first  interim  report  of  the  previous  contract  (Knott  and  Senior,  1973)  ..  It 
is,  in  fact,  the  integral  equation  used  for  H -polarization  in  the  original  computer 

Since  the  primes  now  refer  to  the  point  of  integration,  and  in  view  of  the  present 
definition  of  r,  the  last  terms  on  the  right  hand  sides  of  eqs.  (3.  1)  and  (3.  2) 
have  the  wrong  signs. 
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program  RAMI  A,  and  since  the  kernel  has  at  most  a  first  derivative  singularity, 

the  equation  is  well  suited  to  numerical  solution. 

However,  (2.  9)  is  not  the  omy  integral  equation  that  can  be  obtained.  If, 

for  example,  the  normal  derivative  of  the  field  H  (p)  of  eq.  (2.7)  is  taken  prior 

z 

to  the  observation  point  reaching  C,  we  find 


Y  E1  (s)  -  rj(s)K  <s)  + 
Os  s 


lim  _i_  _o_ 
£  — >  C  4k  On 


i  p(s')ri 

i 

1. 


(1) 

0 


(kr) 


-i(n'  •  r)H^(kr)>  K  (s')d(ks')  .  (2.10) 

1  J  s 

For  the  first  term  in  the  integrand,  the  differentiation  and  integration  processes 
can  be  interchanged  and  the  non-integrable  singularity  that  results  when  £  is  on 
C  can  be  treated  analytically  as  before.  This  is  not  true  of  the  second  term,  how¬ 
ever.  Hence 


YE1  (s) 
0  s 


a  n(3>Ks(S>  -  a 


rj(s’)  KMs'Hn 


r)H^(kr)d(ks') 


lim  _L  JL 

£— *  C  4k  On 


K  (s')  'n'  •  r)H^(kr)d(ks') 

\  s> 

Jc 


(2.11) 


which  is  not  so  convenient  for  numerical  solution  because  of  the  "second 
derivative  singularity"  which  it  contains.; 

We  now  turn  to  the  case  of  E -polarization.  A  particular  advantage  of  the 
impedance  boundary  condition  is  that  it  preserves  the  duality  inherent  in  Maxwell's 
equations  and,  in  the  present  instance,  this  enables  us  to  deduce  the  integral 
equations  for  E-polarization  from  those  we  have  already  given. 
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As  pan  of  the  duality  transformation,  the  surface  electric  current 


K  =  n^H  is  replaced  by  the  surface  magnetic  current 


K  =  -n*E  . 

From  the  boundary  condition  (2,  1),  however. 


(2, 12) 


(2. 13) 


and  similarly,  from  (2,  2), 


K  =  -  Ynn,  K" 
-  !)  0 


(2.  14) 


Hence,  as  regards  the  component  K 


K-- >K 

s  s  0  z 


(2. 15) 


where,  from  (2.3),  K  *  -H  'S  simply  the  axial  component  of  the  surface  elec- 

z  s 

trie  current.  On  applying  the  transformation 

<"*<■  E1^V  Ks  ~ ^  -”Z0K-  ■  Y0-^Z0' 

to  eqs.  \2.  9)  and  (2,  11),  we  immediately  obtain 


Ci)  *  ~  t,(s)K  (s)  +  -  \  K  fs,)Hf(i\kr)d(ks») 
0  z  2  z  4  \  (i 


r 

\ 


~  \  \  r7(s’)Kz(s,)(n' •  r)II^(kr)d{ks’)  ,  (2,16) 

Jc 


H^s)  =  '■^•K  (s)  +  j  X  K  (s*)(n  •  r)H^(kr)d(ks') 
s  2  z  4  z  1 

JC 


£’™C  4k  "k  \  r5(sf)Kz(s’)("'  •  rm^tkrjdflcs*)  .  (2.  17) 


Like  (2,  1)  and  (2. 11),  these  are  alternative  equations  appropriate,  in 
this  case,  to  E-po!arization,  The  first  is  identical  to  that  derived  by  Andreasen 
(1965)  for  incidence  perpendicular  to  the  z  axis,  and  was  quoted  by  Knott  and 
Senior  (1973a;  eq.  3..  2).  Since  the  highest  singularity  of  the  kernel  is  a  first 
derivative  one,  the  equation  is  weli  suited  to  numerical  solution  and  is  the  obvious 
one  to  cnoose  for  E -polarization  In  contrast,  Eq,  (‘l.  17)  has  (in  effect)  a  second 
derivative  singularity  in  the  kernel  and  is  much  less  convenient  to  use,  but  in 
spite  of  this  it  is  apparently  the  equation  employed  in  RAM1A.  This  is  not  en¬ 
tirely  unnatural  since  the  computer  program  was  originally  developed  for  skew 
incidence,  but  even  in  this  more  general  case  it  is  still  possible  to  obtain  two 
coupled  integral  equations  having  only  first  derivative  singularities  at  the  expense 
of  introducing  derivatives  of  the  unknown  currents.  As  Andreasen  (1965)  shows, 
these  latter  derivatives  can  be  eliminated  by  partial  integration,  and  the  resulting 
equations  then  degenerate  to  (2.  9)  and  (2.  16)  for  incidence  pe-pendieular  to  the 
z  axis. 

The  fact  that  the  computer  program  RAMI  A  is  based  on  eqs.  (2.  9)  and 
(2.  17)  explains  our  differing  experiences  in  running  it  for  H-  and  E -polarizations. 
It  is  now  obvious  how  we  can  overcome  the  difficulties  which  formerly  bedevilled 
us  for  E-polarization  and  this  matter  is  taken  up  in  Chapter  V. 

2.  3  Vector  Formulation 

As  we  have  seen,  the  problem  of  a  plane  wave  incident  on  the  cylindrical 
structure  in  a  plane  perpendicular  to  its  axis  can  be  treated  as  a  scalar  one,  but 
since  Hertz  vectors  are  a  convenient  tool  for  any  scattering  problem,  it  is  of 
interest  to  examine  the  integral  equations  resulting  from  this  type  of  vector 
formulation. 

* 

In  terms  of  the  electric  and  magnetic  Hertz  vectors  TT  and  TT  ,  the 
scattered  field  can  be  written  as 


Es  =  vAvArr+ikZ0vAn* 

— S  =  ^L'*-ikY0VAIT 


(2.  18) 
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10 


on  allowing  £  to  approach  the  surface  and  using 

HS(£)  =  K  (s)  -  Hi  (s)  ,  £  on  C  , 

Z  b  2 

the  integral  equation  (2.  9)  is  at  once  obtained.  Similarly  for  the  scattered 

s 

electric  field:  if  the  component  E  is  evaluated  at  a.  point  on  the  surface,  Ibe 

s 

fact  that 

ES(f>)  =  rj(s)Z  K  (s)  -  E*  (s)  ,  £  on  C  , 

leads  immediately  ro  the  integral  equation  (2.  11). 

Thus,  the  vector  formulation  leads  to  precisely  the  same  integral  equations 
as  the  scalar  one.  For  H-polarization,  the  equation  with  the  first  derivative 
singularity  stems  from  a  consideration  of  the  magnetic  field,  whilst  the  electric 
field  gives  rise  to  the  equation  having  the  second  derivative  singularity.  Either 
is  sufficient  to  determine  the  single  non-zero  component  of  the  surface  elec¬ 
tric  current.  If,  instead,  we  had  considered  F~  polarization,  the  electric  field 
would  have  generated  eq.  (2.  16)  and  the  magnetic  field  (2..  17).  This  is  obvious 
from  duality  and  here  again  either  equation  is  sufficient. 

2 ..4  Thin  Sheets 

The  integral  equations  derived  in  the  previous  sections  are  valid  for  a 
two-dimensional  surface  at  which  an  impedance  boundary  condition  is  imposed, 
and  provided  its  boundary  profile  C  encloses  a  region  of  non-zero  area,  either 
of  the  two  equations  (2.  9)  and  (2.  11)  obtained  for  H-polarization  is  sufficient  to 
specify  the  surface  electric  current  K.  This  is  also  true  of  the  E-poIarized 
integral  equations  (2.  16)  and  (2.  17),  and  for  each  polarization,  the  equation  with  the 
first  derivative  singularity  is  the  more  attractive  one  for  numerical  solution. 
However,  as  the  thickness  oi  tne  body  decreases  in  comparison  with  its  length, 
so  does  the  rate  of  convergence  of  a  numerical  solution,  and  as  the  two  sides  of  the 
body  approach  one  another,  these  difficulties  become  severe.  In  fact,  in  the  limit 
of  an  infinitesimally  thin  body,  neither  type  of  integral  equation  is  capable  of 

1 1 


solution  as  it  stands.  This  limit  is  the  case  of  a  thin  sheet  (not  necessarily 
planar),  and  as  we  shall  show,  both  integral  equations  are  now  needed  to  determine 
the  surface  current. 

It  is  again  sufficient  to  restrict  our  attention  to  H-polarization.;  If  the 
r'ylinder  is  infinitely  thin,  the  closed  contour  C  can  be  broken  up  into  two  open 
contours  C+  and  C  lying  on  the  upper  (positive)  and  lower  (negative)  sides,  and 
the  integrals  expressed  as  integrals  over  C  (say)  alone.  Thus,  for  a  thin  sheet 

T 

whose  surface  impedances  are  the  same  on  both  sides,  eq.  (2.  7)  can  be  written  as 


r 


Hz(£)  =Hz(£)+^-  \  •|j2(s')(n^-r)H^(kr)+in(s,)J1(s,)H^1,(kr)id(ks') 

JC  ^  (2.24) 


where 

J  (s)  =  K  (s)+K  (s)  (2.25) 

J  ■  "S  -s 

is  the  sum  of  the  surface  currents  on  the  two  sides  and 

J„(s)  =  K  (s)  -  K  (s)  (2.26) 

*s  -s 

is  the  difference.  J^(s)  and  J2(s)  are  proportional  to  the  local  strengths  of  the 
equivalent  magnetic  and  electric  current  sheets,  respectively.  We  now  take  the 
limits  of  eq.  (2.24)  as  £  approaches  the  surface  successively  from  above  and 
below.  Addition  of  the  two  limits  gives 


Hl(s)  =  iJi(s)+i 


c 


(l) 

0 


(kr  )d(ks’) 


(2.  27) 


which  is  a  simple  integral  equation  for  J^(s).  Subtraction,  on  the  other  hand, 
yields 


T  .  x  i  f  lim  lim 
J2(s)*  J\£-»C  -£-»C 


\  1 

)  \ 

J  Jc 


K 


r )  J  0  (s ')  H  [  ^  (k  r  )d  (ks ') 

&  i 


(2.28) 
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and  by  analytical  examination  of  rhese  limits  ic  can  be  shown  that  (2.  28)  is  merely 
the  identity  ^(s)  =  ^2^'  Unless  ^(s)  =  ^  (and  this  is  not  in  general  true  even 
for  a  uniform  planar  sheet),  the  integral  equation  (2.  9)  has  now  served  to  deter¬ 
mine  only  the  sum  of  the  surface  currents. 

To  find  ^(s)  and.  hence,  the  currents  themselves,  it  is  necessary  to  use 
the  integral  equation  (2,  11)  in  addition  to  (2.  9).  Paralleling  the  analysis  leading 
up  to  (2,  11),  the  normal  derivative  of  (2.24)  is  taken  as  the  observation  point 
approaches  the  sheet  from  above  and  below.  Subtraction  of  the  two  limits  again 
yields  an  identity,  this  time  for  J^(s),  but  addition  produces  the  following  integral 
equation  for  J^s): 

Y0E«(s)  '  2r'{s)32{s)+ 4k  aT  \  Vs’KSyfm'.Vldte*)  .  (2.29) 

+  Jc+ 


This  is  analogous  to  (2.  11)  in  having  a  second  derivative  singularity. 

For  a  thin  sheet  it  is  therefore  necessary  to  employ  both  types  of  integral 
equation  developed  in  Sections  2.2  and  2.3  and  for  H -polarization  the  appro¬ 
priate  equations  are  most  conveniently  expressed  in  the  forms  (2.27)  and  (2.29). 
The  equations  for  E-polarization  are  similar  and  can  be  obtained  using  duality. 
Thus,  from  (2.  27)  we  have 


-E1  (s) 
z 


"  2  J1(S)+  4  \  db) 


(2.30) 


where 


On  writing 


J, (s)  =  -  (E  +E~)  =  rj(s)Z  (H  -H  ) 
1  z  z  0  +s  -s 


J,(s)  =  -q(s)z„ j0(s)  , 


(2.31) 


eq.  (2.  30)  becomes 


Y0Ez(s)  =  2  ^S)J3(S)+  4  \  J3(s’)Hj)1)(kr)d(ks') 

Jc. 


(2.32) 


which  is  an  integral  equation  for  the  strength  of  the  equivalent  electric  current 
sheet.  Similarly,  from  eq.  (2.29), 


V>> 


r 

SM  J4<S)-  £^c+  5  ^  \ 

JC+  (2.33) 


where 


J4(s)  -  (Ez-E2)  =  -J2(s) 


(2.  34) 


is  the  strength  of  the  equivalent  magnetic  current  sheet..  Equations  (2.  32)  and 
(2.34)  are  the  required  integral  equations  for  E-polarization. 

A  special  instance  of  a  thin  sheet  is  a  planar  strip  (or  ribbon),  and  for 
this  geometry  it  is  of  interest  to  pursue  our  study  of  E-polarization  a  little  further. 
If  the  strip  occupies  the  portion  0  x  <  L  of  the  plane  y  =  0,  where  x,y,  z  are 
Cartesian  coordinates,  and  is  illuminated  by  a  plane  wave  having 


i_A  -ik(xcosa+y  sino) 

&  Z  6 


(2.  35) 


eq.  (2.  32)  becomes 


^-ikxcosa^  ir?(x)j^(x)+  I  \  (2.36) 


while  (2.  33)  reduces  to 
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For  a  plane  wave  at  grazing  incidence,  sir  •->  ~  0,  and  eq.  (2.  37)  gives  J^(x)  =  0.. 

We  are  now  left  with  the  single  equation  (2.  36)  t.c  solve,  and  since  E  =  E  ,  its 

z  z 

solution  specifies  the  surface  fields  on  the  strip.  In  this  case  the  impedance  strip 
is  equivalent  to  a  resistive  one  (see  Chapter  IV). 

If  L  =  oo  and  rj  =  constant,  the  strip  is  a  uniform  "metallic"  half  plane. 
The  integral  equations  (2,36)  and  (2,  37)  are  then  identical  to  cqs,  (11)  and  (12) 
of  Senior  (1952)  and  can  be  solved  analytically  using  the  Wiener--Kopf  technique. 
From  the  resulting  expression  for  the  scattered  field,  the  edge  diffraction  co¬ 
efficient  can  be  determined.  This  is  shown  in  the  next  Chapter. 
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EDGE  DIFFRACTION  ,jY  AN  IMPEDANCF  HALF  PLANE 

The  diffraction  of  a  plane  \va  je  by  a  half  plane  at  which  an  impedance 
boundary  condition  is  imposed  has  been  soiveci  for  both  normal  and  oblique  0,  /) 

incidence  (Senior,  1952,  1960),  and  even  for  the  general  problem  of  an  impedance 
wedge  of  arbitrary  included  angle,  solutions  are  available  (Senior,  1959; 
Maliuzhinets,  1960;  Lebedev  and  Skal’skaya,  1963)  for  normal  incidence.  Never¬ 
theless,  computed  data  showing  the  effect  of  the  impedance  on  the  scattering  are 
very  scarce.  Senior  (1952)  oresented  some  patterns  for  a  specific  complex  value 
of  rj  and  Maliuzhinets  (1960)  has  reported  that  one  of  the  key  functions  involved 
in  the  solutions  has  been  tabulated,  but  no  reference  is  given. 

The  analytical  solutions  can  be  used  to  find  the  diffraction  coefficients 
associated  with  the  edge.  Maliuzhinets'  expressions  are  particularly  convenient 
fo"  this  purpose,  and  from  the  resulting  coefficients  for  an  impedance  half  plane. 
Bowman  (1967)  obtained  a  high  frequency  approximation  to  the  backscattered  field 
of  an  absorbing  strip  at  normal  incidence.  The  predicted  values  are  in  good 
agreement  with  measured  data. 

If  the  half  plane  is  illuminated  by  an  E-  or  H-polarized  plane  wave  of  unit 
amplitude  incident  at  an  angle  a  in  a  plane  perpendicular  to  the  edge,  the  cor¬ 
responding  diffracted  field  at  large  distances  from  the  edge  and  viewed  in  a  direc¬ 
tion  6  away  from  either  of  the  optical  boundaries  is 


nr  , 

im e 


where  (Maliuzhinets,  i960) 


a 


sin 


r 


t  (-e) 


(2n-0) 


.  Q  a 

sin  —  +  cos  — 


sin  -  -  cos  ~ 


(3.1) 
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with 


(3.2) 


C (/3)  •(/,  (fi+n+xH  $+n-x)<i  (fi-n+xh  {(3-tt-x)  , 

n  7i  ft  ft 


n  sinv-2i/i~77sin  tj-  +  2v 

^  \  ^  ~ 

(  Jo 


(3.3) 


The  angle  x  is  ^uch  thf 


cc  ?x  = 


l/^  for  E -polarization 
n  for  H-polarization 


'3.4) 


Many  properties  and  alten  itive  definitions  of  the  function  ^  (/3)  are  listed  by 
Bowman  (1967)  in  the  Appet  dix  to  his  paper. 

The  edge  diffraction  :oefficient  P(a,  9)  is  a  function  of  tne  angle  of 
incidence  a,  and  the  case  of  edge-on  incidence  ( a  =  n)  is  of  particular  interest 
to  ua.  We  first  consider  bad  scattering,  d  =  as7r,  and  then  examine  the 
more  general  case  of  bistatic  ;  cattering. 

3.  i  Backscattering  for  Edge-or  Incidence 

From  eq.  (3. 1)  with  9  =  a  -  it ,  we  have 

m  <3-5> 

since  0(0)  is  an  even  function.  As  shown  by  Bowman  (1967), 


^0/2)]“  ,  „  ^ 


so  that 
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and  since 


it  follows  that 


, ,  i  K <ir/2)}8 

b(ir)  =  gCOSX  - - - ; 

c(0)  =  {^Or+xH  ^(tt-x)  j-2  , 


P(n,n)  =  -  y^rcosx 


K<*/2>)f 


{iijr(jr+x)^{r-x)l 


But 


'  2  M 


2  VC0S  4  +  COS  2 


and  hence 


P(7r,  7r)  =  -i 


4cosx 


V  (x) 

71 


1+^cosf  VlV^2’} 


(3.6) 


these. 


Bowman  has  derived  two  alternative  forms  for  c  (/3).  From  the  first  of 

7 r 


V**  -1/8  /  --  „V/2 

Mir/2)  *  <8cosx)  I1  +1/2  COS  exp 


x'' 


ICO 


vdv 
in  \  cos  v 


J  (3. 7) 


giving 
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C  ^ioo 


n/  '  -  1  I  2  \  V  dv  ( 

P(ff,  vi  -  -  -  exp<  -  \ - > 

-  7T  \  COS  V 

L  Jx  y 


(3.8) 


whereas  from  the  second 


rt/2)  -  2“1/2  (l  +  /2  cos  |)  exp  -j 


r  p^/2 

1  '  7rsinv-2v 


8?r 

V.  JX 


COS  V 


giving 


dv 


(3.9) 


P(7r,  n)  =  -  7  cos  x  exp  J  -  - 
4 


r  p?r/2  -x  '"N 

7TCOS  U  -  fa  -  2u) 
sin  u 

'0 


du 


(3.10) 


Equations  (3.  8)  and  (3. 10)  can  be  used  to  obtain  analytical  approximations 
to  P(ir,Tr)  for  small  and  large  |rj  j..  For  this  purpose,  consider  E-polarization  tor 
which  cosx  -  1  /n.  If  |r?|  «  1, 


_1 

n 


implying 


ewise 


X  -  -i  log 


R 

2 


v 

cos  V 


2ve 


iv 


and  therefore 


(3.11) 


/">  l  OD 


vdv 

COS  V 


--2(l-iX)eiX  - 


n(i  -  log  ^ ) 


from  which  we  obtain 
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P(jr,  7r)  - 


When  n  -  0,  eq.  (3. 12)  gives 


P(?rJ7r)  =  -  | 


(3.12) 


(3.13) 


in  agreement  with  the  known  result  for  a  perfectly  conducting  half  plane.  If,  on  the 
other  hand,  |r)  j»  1.  we  have 


X  - 


v 


(3. 14) 


and 


so  that 


-X 

7TCQS  u-  (?r-2u) 
sinu 


i  /  2  \ 

6XPC 


(3.  15) 


To  see  how  successful  the  approximations  (3. 12)  and  (3. 15)  are  in  covering 
the  entire  range  0  v  |p  K  oo,  P(n,  tt)  has  been  computed  for  real  r?  by  direct 
numerical  evaluation  of  one  of  its  integral  expressions.  Since 


uoo 


vdv 
cos  v 


=  -2k 


where  k  =  0.  9159656. .. ..  is  Catalan's  constant  (Bowman,  1967),  eq.  (3  8)  can  be 
written  as 


P(?r,  7r)  =  -  ~  exp  -j  -  ~ 


(3.16) 
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This  integral  was  evaluated  using  a  32  point  Gauss  quadrature  formula  and 
Fig.  3-1  shows  a  plot  of  the  normalized  diffraction  coefficient 


P(tt,  7 r) 


P  (tt  j  7f) 


Pfcr,  n) 


n=o 


(3.17) 


for  real  n,  where  P(7r,  tt)  j  is  given  in  eq.  (3.13).  Results  computed  using 

>ri-  0 

the  approximate  formulae  are  included  for  comparison  and  it  will  be  observed  that 
(3. 12)  and  (3. 15)  cover  the  entire  range  of  |r?  |  with  surprising  accuracy:  even  for 
|r?j  =  1  the  large  and  small  impedance  approximations  differ  by  only  20  percent  and 
are  individually  in  error  by  no  more  than  about  10  percent.  Equations  (3. 12)  and 
(3.,  15)  are  therefore  adequate  for  all  practical  purposes  and  these  are,  of  course, 
not  restricted  to  real  values  of  i). 

The  corresponding  result"  for  H-polarization  follow  on  replacing  n  by  l/rj. 


3. 2  Bistatic  Scattering  for  Edge-on  Incidence 

The  analysis  in  this  more  general  case  is  not  dissimilar  to  the  one  we  have 
already  carried  out,  although  the  expressions  obtained  are  more  complicated. 
From  eq.  (3. 1)  with  a  -  tt,  the  bistatic  diffraction  coefficient  is 


8 

.  cosec  -- 

POr,  0)  =  -  7 


U<0) 


^c-(-0)+V'(2jr-0)  J- 


(3.18) 


which  can  be  written  as 


P(7T,  8)  =  -  j-  cosec  | 

2  2  cosx+sm$  i/(0) 


(3. 19) 


on  using  the  relations  (A.  14)  and  (A.  15)  of  Bowman  (1967).  Using  relation  (A,  13) 
of  the  same  reference,  it  can  be  shown  that 


m 


cos 


l+uosg+fj  f  1  2 


X-,  / tt  0\  «  u  (ir+x)(t  (tt-x)  ! 

COS  £  +  COS(j  -  -J  TT  A  TT  J 


4i 


|  +  CO»(f  +  f)  / 


C0SX  C03  |  +  COS  1^-  -  -y 


f,  6  \  V00S  2  +  cos 


r  ~\  2 

_v>  (8-tt+xH  (e-TT-x) 

1  Pfa.ri 


<»  (x></  (-x) 

7T  7 T 
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and  hence 


The  exponent  is  a  finite  range,  real -variable  integral  which  converges  for  all  x  and 
6,  and  can  be  evaluated  nunericaliy.  The  right  hand  side  of  eq.  (3.21)  is  unity  for 
6=0,  as  required,  but  tends  to  infinity  as  6  —*ir  regardless  of  x-  This  is  the 
forward  scatter  case  in  which  the  direction  of  observation  coincides  with  the  reflected 
wave  and  shadow  boundaries,  and  the  infinity  is  attributable  to  this  fact. 
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Equation  (3.21)  can  be  used  to  obtain  analytical  approximations  to 
P(7 r.ir-b)  for  small  and  large  Jp  {.  If  Jr? ) «  1,  x  is  given  approximately  by  (3. 11) 
and 


r  r  6  ^ 

P  (7T.7T-6)  6_  I  1  \  iX 

P(tt,  71)  '  sec  2  exp)  2^  \  {_i)(  e  cos  u)du  \ 

^  Jo  J 

6 

-  sec  —  exp 

Hence,  from  (3. 12), 

P(ff,  n-  6)  —  -  ~  sec  ~  exp  -i  -  ^  j^2  -  (1  +cos  6)  log  ^  j  ■  (3.22) 

Not  surprisingly,  this  is  consistent  with  eq.  (3.  12)  when  6=0,  and  it  also  reduces 
to  the  known  bistatic  diffraction  coefficient  for  a  perfectly  conducting  half  plane 
when  n  -  0.  For  real  p,  P(tt,  tt-6)  increases  monotonically  as  6  increases  from 
0  to  ; r. 

If  |r?|»  1,  X  can  be  approximated  by  the  value  given  in  eq.  (3. 14),  but 
in  order  to  arrive  at  an  expression  for  the  diffraction  coefficient  it  is  necessary 
that 


-  ^  (1  -  cos  6)  log  ^ 


1 


(3.23) 


With  this  restriction. 


P(7r,7r-6)  ^  6 

P<7T,  7 r) 
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^  r  6 
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6 

sec  - 


j  1  /  6  ^  j 

-  exp  <  1--—  j  i 

4  I  Trp  \v  sin  6/  [ 
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and  on  inserting  the  expression  for  P(7r,  7r)  given  in  eq.  (3. 15),  we  have 
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.w  .  ^  1  6  I  1  /  6 

P(tt,  7T-6) - —  sec  ~  exp  •<  -  —  11  +  - r 

■in  2  7rrj  \  sin  6 


(3.24) 


This  is  in  agreement  with  eq.  (3. 15)  when  6=0  and,  consistent  with  (3.23), 
remains  small  as  6  increases. 

The  small  and  large  impedance  approximations  (3.22)  and  (3.24)  both 
indicate  that  P(r, n-6)  changes  rather  slowly  as  6  increases  from  zero,  and  this 
has  been  confirmed  by  a  numerical  determination  of  P(7r,  r-  6)  using  an  exact 
integral  expression  analogous  to  that  in  eq.  (3.  21). 


IV 


INTEGRAL  EQUATIONS  FOR  THIN  SHEETS 

One  of  the  main  advantages  of  the  impedance  boundary  condition  is  its 
ability  to  simulate  a  variety  of  surface  materials  and/or  imperfections  and  to  do 
so  in  terms  of  a  single  parameter  only.  Unfortunately  this  generality  is  also  a 
weakness,  and  having  used  the  integral  equations  or  some  other  theoretical 
approach  to  arrive  at  a  desirable  specification  of  the  surface  impedance,  it  is 
by  no  means  evident  how  (or  even  if)  this  impedance  can  be  realized  in  practice. 

To  avoid  the  difficulty,  it  is  natural  to  seek  a  formulation  in  which  the 
properties  of  any  surface  or  coating  are  included  explicitly,  both  geometrically  and 
electrically..  Such  a  formulation  has  been  considered  by  Osniro  et  al.  (1971),  who 
developed  the  integral  equations  for  the  problem  of  an  E-polarized  plane  wave  in¬ 
cident  on  a  metallic  cylinder  fully  or  partially  clad  with  a  dielectric  layer.  The 
geometry  and  permittivity  of  the  layer  are  arbitrary,  but  its  permeability  is 
assumed  to  be  that  of  free  space.  Since  the  field  scattered  by  the  layer  is 
expressed  as  a  surface  integral  of  the  polarization  currents,  the  resulting 
integral  equations  are  difficult  to  solve  in  general,  but  are  amenable  to  solution 
in  the  limiting  case  of  an  infinitesimally  thin  layer.  The  computer  program  based 
on  this  formulation  is  designated  REST  and  has  proved  very  helpful  in  our  studies. 

Although  the  problem  of  a  plane  wave  incident  on  a  resistive  sheet  in 
isolation,  i.  e. ,  without  the  scatterer  present,  is  of  little  practical  interest  in 
itself.  It  is  instructive  to  compare  it  with  the  problem  analyzed  in  Section  2.4. 

For  this  purpose,  attention  is  confined  to  E-polarization,  for  which  program  REST 
is  applicable. 

4.  1  Electrically  Resistive  Sheets 

Consider  a  thin  sheet  of  highly  conducting  material  whose  permeability 
p  is  that  of  free  space.  If  cr  is  the  conductivity  and  A  is  the  thickness,  we  can 
define  a  surface  resistivity  R  as 
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1 


(ohms) 


oA  ue_x  A  '  •  11 

0  e 

where  x  is  the  electric  susceptibility  and  en  is  the  permittivity  of  free  space, 
e  U 

and  as  A — >0  we  can  imagine  cr  to  be  increased  in  such  a  manner  that  R  is 
finite  and  non-zero  in  the  limit.  The  result  is  an  idealized  (infinitesimally  thin) 
electric  resistive  sheet  whose  electromagnetic  properties  are  specified  by  the 
single  measurable  quantity  R. 

Mathematically  at  least,  the  sheet  is  simply  an  electric  current  sheet 
whose  strength  is  related  to  the  tangential  electric  field  via  the  resistivity.  Since 
p  =  /Uq,  there  is  no  magnetic  current  present  and 

V(E+-E)=0  (4,2) 

where  the  affices  t  refer  to  the  upper  (positive)  and  lower  (negative)  sides  of  the 
sheet  and  n  =  n+  is  the  outward  normal  to  the  upper  surface.  If  J  is  the  total 
electric  current  flowing  in  the  sheet,  i.e. ,  its  strength, 

Vh+~h”)=J  (4.3) 

and  from  the  definition  of  the  surface  resistivity, 

+ 

nA  (nAE~)  =  -RJ  ,  (4.  4) 

With  R  specified  (it  may  be  a  function  of  position  on  the  sheet),  eqs.  (4.  2)  through 

(4.  4)  constitute  the  boundary  conditions  at  the  surface. 

The  field  occurring  in  each  of  the  above  equations  is  the  total  (incident  plus 

X  j 

scattered)  field,  and  since  H_  =  H_  when  J  =  0,  eq..  (4.3)  implies  for  a  planar  sheet 


at  least. 


nA  if  = 


Hence,  from  h. 4),  the  surface  current  density  on  the  upper  surface  is 

H+'"a{h‘-8a(s-i)) 


27 


and  If  a  surface  impedance  q  is  attributable  to  the  sheet,  eq,  (4.  5)  suggests  that 

q  =  2YqR  .  (4.6) 

Nevertheless,  the  problems  of  a  resistive  sheet  and  a  sheet  at  which  an 
impedance  boundary  condition  is  enforced  are  different  from  one  another  and,  in 
general,  not  equivalent.  Although  there  is  some  similarity  between  the  conditions 
imposed,  an  impedance  boundary  condition  specifies  a  connection  between  the 
tangential  components  of  .E  and  H  on  each  side  of  the  sheet  and  forces  a  mag¬ 
netic  current  to  occur  in  addition  to  an  electric  one.  The  sheet  is  therefore  a 
combined  electric  and  magnetic  current  sheet  with  the  two  strengths  related  via 
the  impedance.  In  contrast,  a  resistive  sheet  supports  no  magnetic  current  by 
virtue  of  eq.  (4.  2),  Because  of  this  the  two  problems  are  basically  distinct, 
and  unless  the  geometry  is  such  that  no  magnetic  current  is  excited  on  the 
impedance  boundary  condition  sheet,  it  is  not  possible  to  relate  R  and  q. 

To  obtain  the  integral  equations  appropriate  to  a  resistive  sheet,  we 

consider  first  the  case  of  E -polarization.  If  E1  =  zE1,  then  J  =  z  J  and  the 

—  z  —  z 

boundary  condition  (4.4)  becomes 

4- 

E"  =  R  J  (4.7) 

z  z 

where  R  is  the  surface  resistivity,  assumed  known.  It  is  a  trivial  matter  to 
derive  an  integral  equation  for  the  strength  of  the  equivalent  electric  current 
sheet,  and  the  vector  formulation  is  (perhaps)  even  more  straightforward  than  the 
scalar  one.  From  the  definitions  (2.  19)  and  (2.  20)  of  the  electric  and  magnetic 
Hertz  vectors  describing  the  scattered  field,  we  have 

TT  '  (£)  =0  {*.  8) 

on  using  (4.2)  .  Also 

TT(£)  -  -  ~~  z  \  Jz(s')H^(kr)ds’  (4.9) 
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where  the  integration  is  along  the  sheet  (or  on  one  side  only),  and  hence,  from 
(2.  18), 


Z, 


E,(£>  =  e!,(£)  -  -J-  \  Jz(s,)H^1)(kr)d(ks'). 


(4.  10) 


If  £  is  now  chosen  to  lie  on  the  sheet,  application  of  the  boundary  condition  (4.7) 
gives 


Y0Ez(s)  =  Y0R(s)Jz(s)  +  i  \  Jz(s’)HQl)(kr)d(ks')  ,  (4.11) 


which  is  a  rather  simple  integration  equation  for  the  total  current  J  (s).  Note 

z 

that  consideration  of  the  magnetic  field  would  have  produced  only  an  identity  for 
Jz(s).  Alternatively,  the  scalar  formulation  can  be  used  starting  with  eq.  (2.  6) 
and  taking  u  to  be  the  z  component  of  the  scattered  electric  field.  The  resulting 
integral  equation  is  identical  to  (4. 11). 

From  eq.  (2.31), 

J  (s)  =  - (H  -H  )=J,(s)  (4.12) 

z  +s  -s  3 

where  J^s)  is  the  total  current  satisfying  the  integral  equation  (2.32).  Com¬ 
parison  of  (4.  11)  and  (2.  32)  shows  them  to  be  identical  if  R  and  r?  are  related 
through  eq.  (4.  6)  With  this  identification,  the  electric  currents  are  the  same 
for  resistive  and  impedance  boundary  condition  sheets,  but  since  the  latter  sheet 
can  also  support  a  magnetic  current,  the  scattered  fields  will  differ  in  the  two 
cases  unless  the  magnetic  current  is  not  excited.  A  situation  where  this  occurs 
was  noted  in  Section  2.4  and  hence,  for  a  planar  sheet  at  edge-on  incidence, 
a  resistive  sheet  is  entirely  equivalent  to  an  impedance  boundary  condition  sheet 
whose  surface  impedance  is  given  by  eq.  (4.  6)  . 

« 
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This  conclusion  is  important  in  enabling  us  to  use  the  analytical  information 
available  about  vho  impedance  boundary  condition  sheet  in  studying  the  properties 
of  resistive  sheets.  From  the  exact  analytical  solution  for  a  constant  impedance 
half  plane,  the  edge  diffraction  coefficent  can  be  obtained  as  shown  in  the  previous 
chapter.  For  the  pairicular  case  of  edge~on  incidence,  this  coefficient  is  also 
applicable  to  a  resistive  sheet  for  all  angles  of  diffraction  provided  r?  =  2Y^R 
where  R  is  the  suiface  resistivity  in  a  neighborhood  of  the  edge. 

Turning  now  to  H-pularization  for  which  11*  =  zh\  the  strength  of  the 
resulting  current  sneet  is  J  =  s  and  the  boundary  condition  (4,4)  yields 

+ 

E  =  RJ  .:  (4.  13) 

s  s 

For  the  field  scattered  by  the  sheet  the  magnetic  Hertz  vector  TT  (£)  is  again  zero, 
and  since 

z  r 

rj(£)  =  \  s‘Js(s‘)*y!  Grids'  ,  (4.14) 

g 

it  is  a  trivial  matter  to  find  E  (£),  When  £  is  on  the  sheet,  however,  eq,  (4, 13) 
implies 

ES  ~  -E1  +R  J 
s  s  s 


and  hence 


V>  =  yor(s)‘Vs)+A 


_l_  JL 

4k  9n 


J  (s')(n'-r)H 
s  1 


(kr)d(ks’) 


(4.  15) 


This  is  an  integral  equation  from  which  to  determine  J  (s).  Unfortunately,  it  is 

s 

characterized  by  the  same  type  of  second  derivative  singularity  that  we  have  met 

oefore  and  which  is  difficult  to  handle  numerically,  but  in  suite  of  this,  eq.  (4. 15) 

is  the  one  that  must  be  solved.  Indeed,  tne  consideration  of  the  magnetic  field 

scattered  by  the  sheet  produces  only  an  identity  for  J  (s) . 

s 
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Even  for  a  resistive  sheet  we  are  now  faced  with  two  integral  equations 
analogous  to  the  ones  obtained  for  an  impedance  boundary  condition  sheet,  but 
whereas  the  latter  requires  both  equations  to  find  the  two  currents  that  either  an 
E-  or  H-polarized  wave  excites,  the  resistive  shee:  needs  only  a  single  equation 
for  each  polarization.  Though  one  is  convenient  for  numerical  solution,  the  other 
(for  H-polarization)  is  not.  Were  we  to  compare  eq.  (4. 15)  with  the  analogous 
equation  for  an  impedance  boundary  condition  sheet,  we  would  again  find  them  to 
be  identical  if  17  =  2YqR. 

4.  2  Generalized  Resistive  Sheets 

Our  attention  thus  far  has  been  concerned  w'ith  the  concept  of  an  electric 
resistive  sheet  whose  electromagnetic  properties  are  specified  entirely  by  the 
surface  (electric)  resistivity  R,  but  we  can  also  conceive  of  the  electromagnetic 
dual  of  such  a  sheet.  This  could  be  simulated  using  a  material  having  c  = 
and  high  magnetic  loss,  and  we  could  then  define  a  surface  "magnetic  resistivity" 


R 


i 


(mhos) 


(4.  16) 


analogously  to  eq.  (4.  1)  ,  where  x^  is  the  magnetic  susceptibility,  is  the 
permeability  of  free  space  and  A  is  the  thickness  of  the  layer.  If  xm  — **00  (or, 
more  realistically,  ,jlraxj  as  A — >0  in  such  a  way  that  the  product  remains 
finite,  we  are  then  left  with  an  idealized  (infinitesimally  thin)  magnetic  resistive 
sheet  whose  properties  are  completely  specified  by  the  single  measurable  quantity 
R  . 

For  a  magnetic  resistive  sheet,  the  boundary  conditions  are  simply  the 
dual  of  those  in  the  previous  section  and  are 

nA(H+-H_)  =  0  (4.  17) 

n  A  (E  +  -  E~)  =  -j"  (4.18) 

nA(nAH')  =  -rV  (4-19> 
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where  the  affices  t  again  refer  to  the  upper  (positive)  and  lower  (negative)  sides 
of  the  sheet  and,  say,  n  =  n  .  The  sheet  supports  only  a  magnetic  current  J 
and,  as  a  result,  the  tangential  components  of  the  magnetic  field  are  continuous 
across  it. 

Integral  equations  from  which  to  compute  J  are  immediately 

obtainable  from  those  in  Section  4.  1  by  invoking  duality.  Thus,  for  H-polari- 

zation,  H1  =  zH1  and  J  =  zJv.  Hence,  from  eq.  (4.  11), 
z  z 

n 

Z0H^(s)  =  Z0R'‘(s)j*(s)  +  |  \  J^(s (kr)d(ks ')  ..  (4.20) 

Similarly,  for  E-polarization,  E1  =  z  E1 ,  J  =  s  J  '  and,  from  (4.  15), 

—  z  -  s 


Z  H1  (s)  =  Z  R  (s)J*  (s)  + 
Os  0  s 


lim  _l _ 

£— >  C  4k  9n 


(s')(n’-r)H 


(1) 

1 


(kr)d(ks').: 


(4.21) 


These  are  analogous  to  the  equations  for  the  magnetic  current  supported  by  an 
impedance  boundary  condition  sheet  having  surface  impedance 


n  =  y0/(2e  v)  , 


(4.22) 


though  we  remark  that  the  problems  of  impedance  boundary  condition  and  mag¬ 
netic  resistive  sheets  are  completely  equivalent  only  in  those  cases  where  an 
impedance  boundary  condition  sheet  supports  a  single-component  current  alone.; 
One  such  situation  is  a  planar  sheet  at  edge-on  incidence. 

ft  is  now  logical  to  consider  a  generalized  resistive  sheet  consisting  of 
a  superposition  of  an  electric  and  a  magnetic  resistive  one.  With  such  a  sheet, 
the  electric  and  magnetic  currents  are  produced  by  the  electric  and  magnetic 
susceptibilities  of  the  material  respectively,  and  are  specified  by  R  and  R  ' 
independently  of  one  another.  From  eqs.  (4.2)  through  (4.4)  and  (4.  17) 
through  (4. 19),  the  electromagnetic  boundary  conditions  at  a  generalized  resis¬ 
tive  sheet  are 
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(4.23) 

A  ,  +  -  * 

nA(E  -E  )  =  -J 

(4,24) 

iiA{SA(E++£")}  =  -2RJ 

(4,25) 

n  A  |n  A  (H_+ + H  )j  =  -2R  "j  ' 

(4.26) 

where  the  last  two  have  been  deduced  from  (4.  4)  and  (4,  19)  by  writing  these  as 

conditions  on  the  sum  fields.  This  avoids  the  cross-coupling  occurring  at  each 
side  of  the  sheet. 

The  Hertz  vectors  describing  the  scattered  fields  are 


=  “  4k  \  ifcOH^Mds'  (4.27) 


U"(e>  =  - 


4k 


*  (1) 

J  (s»)Hg  '(kr)ds* 


(4.28) 


s  s 

from  which  E_  (£)  and  II  (p)  can  be  computed  using  the  eqs,.  (2.  18).,  By  taking 
the  observation  point  to  lie  on  the  sheet  and  imposing  the  boundary  conditions 
(4.  15)  and  (4.  16),  we  then  obtain  integral  equations  from  which  J  and  _J*  can 
be  computed.  Though  it  is  natural  to  expect  that  the  equations  will  be  coupled, 
this  is  not  in  fact  the  case. 

We  shall  again  confine  attention  to  the  principal  polarizations  and  it  is 
actually  sufficient  to  consider  E -polarization  only.  If  E1  =  zE* ,  then  J  =  z  J 


and  J 


* 


s  J  .  From  the  eqs.  (2.  18),  the  fields  scattered  by  the  sheet  are 
s 


HS<£)  - 

S 


JL 

4k  9n 


J%,)(n'-r)H(11)(kr)d(ks') 
s  1 


4  \  J  (s’)(n  r)H^(kr)d(ks’) 

4  |  z  1 


(4.29) 
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E> 


J  (s'lH^toldlks') 
z  0 


7  \  J  ,!(s’)(n‘ •  r)H^Ctr)d(ks') 

4  v  s  1 

Jc 


(4.30) 


If  the  observation  point  £  ’s  now  chosen  to  lie  successively  on  the  upper  and 
lower  surfaces  of  the  sheet,  the  resulting  fields  can  be  inserted  into  „he  boundary 
conditions 


H++H~  -  2rV'!  , 
s  s  s 


E  4  E  =  2RJ 
z  z  z 


(4.31) 


Since  the  contributions  of  the  second  integrals  in  (4.  29)  and  (4.  30)  cancel,  we 
are  left  with 


Z/W  =  Z  R  ^  (s)j  (s)  +  ^  ~  ~ 

Os  0  s  £-^>C  4k  dn 


J  (s')(n' •  r)H^(kr)d(ks*)  . 
s  1 


(4.  32) 


YnEX(s)  *  Y„R(s)J  (s)+  j  \  J  (s')H^(kr)d(ks') 
0  z  0  z  4  \  z  0 


(4.33) 


which  are  two  uncoupled  integral  equat  ions  for  J  and  J  .  They  are  identical 

s  z 

to  the  equations  (4.  21)  and  (4.  11)  ior  individual  magnetic  and  electric  resistive 
sheets,  and  the  scattered  fields  (4.  29)  and  (4.  30)  are  simply  the  sums  of  the  fields 
scattered  by  the  separate  sheets.;  Though  a  generalized  resistive  sheet  does  bear 
some  resemblance  to  an  impedance  boundary  condition  sheet,  it  should  be  noted 
that  in  the  present  case  J  and  J*  are  specified  independently  by  R  and  Rv 
respectively*  the  constraint  (2. 13)  or  (2.  14)  which  characterizes  an  impedance 
boundary  condition  surface  would  only  exist  if  R  =  1/ (4R  ) . 
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An  important  practical  consequence  of  the  above  results  is  that  for  a 
generalized  resistive  sheet  it  is  sufficient  to  solve  only  the  two  integral  equations 
(4.  11)  and  (4. 15)  for  an  electric  resistive  sheet  with  E-  and  H-polarizations, 
respectively.  The  analogous  equations,  (4.21)  and  (4.22),  for  a  magnetic 
resistive  sheet  follow  by  invoking  duality;  and  for  a  generalized  resistive  sheet, 
tne  integral  equations  are  the  same  as  for  a  superposition  of  the  two  sheets,  as 
are  the  scattered  fields.  Of  course,  we  are  still  faced  with  the  numerical  diffi¬ 
culty  posed  by  the  second  derivative  singularity  in  eq..  (4.  15),  but  before  turning 
to  this,  it  is  of  interest  to  consider  the  problem  of  a  generalized  resistive  sheet 
in  the  presence  of  a  body. 


4.  3  Sheets  and  Bodies 

A  resistive  sheet  is  of  concern  to  us  because  of  its  potential  for  reducing 
the  scattering  from  a  body  in  its  vicinity.  Until  now  our  studies  have  been 
confined  to  an  isolated  sheet,  and  though  these  have  served  to  pinpoint  the  types 
of  integral  equation  which  are  appropriate,  the  practically  important  problem 
is  that  of  a  sheet  in  the  presence  of  a  body.  This  is  the  problem  that  we  now 
address  and,  as  we  shall  see,  the  analysis  is  a  rather  trivial  extension  of  that 
which  has  been  done  before 

Consider  a  body  whose  profile  is  the  (closed)  contour  C ^  and  which  has 
the  impedance  boundary  condition  (2,  1)  or  (2.  2)  imposed  at  its  surface.  Some¬ 
where  in  the  vicinity  of  the  body  there  Is  a  generalized  resistive  sheet  whose 
profile  is  the  open  contour  and  which  is  subject  to  the  boundary  conditions 
(4.23)  through  (4.26).  A  plane  electromagnetic  wave  is  incident  and  it  is  suffi¬ 
cient  to  assume  this  to  be  E-polarized,  i.  e. ,  E1  =  zE1 .  Hence,  E  =  zE  and 

Z  Z 

ll^xH^+yH^,  and  from  the  eqs.  (4.23)  through  (4.26), 


J  =  zj 


J  —  S  J 


(4.  34) 


where  J  and  J  ‘  are  the  electric  and  magnetic  currents,  respectively,  of  the 
equivalent  current  sheet.  The  Hertz  vectors  describing  the  field  scattered  by 
the  sheet  are  therefore 
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U  (£)  - 


A^O 
Z  4k 


Jz(s,)H“)<kr01)ds'  ■ 


(4.35) 


1 


ZJ(1);::(£>  - 


Ik  \  S'Js(sW(0\r01)ds’ 


(4.  36) 


where  r^.  -  £-£!  ,  i  =  1  or  2. 

The  body  supports  surface  currents  K  and  K'v  where  (see  eqs.  2.  3  and 

2.12) 


K  =  zK  .  K  =&K 
z  —  s 


(4.37) 


Since  Kg  -  -rjZ^K^  by  virtue  of  the  boundary  condition,  the  Hertz  vectors 
describing  the  field  scattered  by  the  body  are 


U(2)(£> 


„!o 

Z  4k 


K  (s’)H^(kr  )dsT 
z  u  U2 


(4.  38) 


U(2);:(£) 


\  s  *  rj(s ')  Kz  (s  *) Hq1  }  (kr Q2 ) ds  ’ 

JC„ 


(4.39) 


The  complete  scattered  field  is  given  by  ..tie  eqs.  (2.  18)  with 

U  =  n(1>+n<2>,  TT*=Tr1,SW2)*  , 

and  at  an  arbitrary  point  £  in  space 
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~~T  \  Js(s,)(",^01)H(i1)(kr01)d(ks«) 
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4  \  Kz(s,)H01)(kr02>d(ks,) 
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+  4  \  r7(s»)Kz(s')(n'-?02)H(ll)(kr02)d(ks‘)  , 

Jc~ 


(4.40) 


H(£)  -n\e)  +  (x~  -y^ 


Li  \  JZ(s,)Ho1)(kr01)d(ks*) 
JC1 


0  *, 


,(1), 


4k  \  Js(s,)(n,'r01)Hi  (kr01)d(ks«) 


+  ^  \  Kz(s%^(kro2)d(kS0 


+  4k  \  HfeOK  <8*)(n'-r  )H(1)<kr  JdfltsO 


Jc, 


(4.41) 


We  observe  that  there  are  three  unknown  components  of  the  currents,  two  for  the 
sheet  and  one  on  the  body,  ano  three  integral  equations  are  therefore  necessary 
to  determine  them. 

Let  £  approach  a  point  s^  on  the  sheet  successively  from  above  and 
below.  From  the  second  of  the  eqs.  (4.  31)  and  the  fact  that 
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lim  lim 

£-^C1+  £-»Cr 


Jg(s,)(n*  •  r01>H1  (kr01)d(ksf)  =  0,  (4,42) 


we  have 


iC 


.(1), 


Y0Ez(si)=Y0R(si)Jz(si)+i  \  Jz(s‘)H0i/(kril)d(kS,) 

Jci 


+  7  \  K  (s»)H!.1)(krl9)d(ks*) 
4  ,  z  0  12 

jC2 


4  \  n(s,)Kz(s,)(n*  •r12)H^(kr12)d(ksl) 

)c„ 


(4.43) 


where  r..  =  £.-£!,  i,j  =  1  or  2.  Similarly,  by  considering  the  magnetic  field 
component  s  •  H  and  using  the  first  of  the  eqs.  (4.  31), 


i  .  .  „  *.  . , ❖.  .  ,  lim  0  b 

H  (s.)  =  R  (s  )J  (s  )+  —  7— ■ 

s  i  1  s  1  p— >  £  4k  dn 


J  (s'Mn' -r  )H^(kr  )d(ks’) 
s  01  1  01 


+  1  \  '■•/b')(fi1f12>H(11)(kr12)d(ks-) 


1 


4T  \  r,(s’)KZ{s,)  {(fi,-?l2)H(ll)(kri2)|d(kS,)- 

C2  (4.44) 


Equations  (4.43)  and  (4.44)  are  necessary  for  the  determination  of  J  and  Jg  .; 
The  first  integral  on  the  right  hand  side  of  (4.  44)  has  a  second  derivative  singu¬ 
larity  of  the  type  referred  to  earlier.  The  last  integral  is  similar  in  form,  but 
the  integrand  is  actually  free  of  singularities  if  and  Cg  have  no  point  in 
common,  and  the  differentiation  can  be  carried  out  explicitly. 
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The  third  (and  last)  of  the  integral  equations  required  is  obtained  by 

allowing  £  to  approach  a  point  on  the  body..  However,  there  are  two  forms 

that  this  equation  can  take.  If  we  consider  E  and  use  the  fact  that  E  =  nZ.K 

z  z  0  z 

on  C 2*  we  find 


Y0£'1Z(S2)  =  2  n<S2)Kzfe2)  +  4  \  KZ(S?)H0  {kr22)d(kS') 


r)(s,)Kz(s,)(nf  •  r22)Hi  (kr22)d(ks0 


+  4  \  ^  z (s'JHq1^ (kr2 i)d(ks') 


+  “  \  J^s’Hn’  . 


(4.45) 


Alternatively,  from  the  magnetic  field  component  s  -H, 


Hs(S2}  =  “  2  Kz(S2)  +  4  \  Kz(s?)("2  ?22)H1  (kl22)d(ks') 


hm  *l(s*)K  (s*)(fi'-r  JH^Ocr  Jd(ks») 

£— >  £2  4k  dn2  \  z  02  1  02 

^C2 


+  4  \  Jz(s')(a2-?2l)H(1,(kr21)d(ks') 


+  Ik  \  J>’>  |e' ' f2i)H<i1)(kr2i>[ d<ks,)  '  (4- 46) 
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This  last  equation  is  similar  to  (4.44)  and  the  same  comments  apply.  If 

J  (s ')  =  0  =  J  (s'),  i.e. ,  if  there  is  no  resistive  sheet  present,  eqs.  (4.45)  and 
z  s 

(4.46)  reduce  to  (2. 16)  and  (2. 17),  as  required. 

For  the  complete  problem  of  a  generalized  resistive  sheet  in  the  presence 

of  a  body,  the  integral  equations  needed  are  (4.  43),  (4. 44)  and  one  of  (4. 45),  (4.  46) 

Because  of  its  greater  simplicity,  (4.45)  is  the  natural  one  to  choose.  If  the  sheet 

has  only  an  electrical  resistivity,  so  that  R  '  is  infinite,  implying  J'"  =  0,  and  if, 

s 

in  addition,  the  body  is  perfectly  conducting  (rj  —  0),  eqs.  (4.43)  and  (4.  45)  become 


Y0Ei<S2> 


J  >d(ks’>  +  7 

z  0  Z1  4 


k  (si)h; 

z  0 


(kr12)d(ks*)  , 
(4.48) 


and  these  are  the  ones  used  in  program  REST.  However,  one  of  the  great 
advantages  of  the  impedance  boundary  condition  is  its  preservation  of  duality, 
which  allows  us  to  deduce  the  results  for  H-polarization  by  applying  the  follow¬ 
ing  transformation  to  the  E -polarized  equations: 


K  - ^  —  rjZ  K  ,  R  4 — >  R  ,  r?4— »l/n,  Yn  Zn  . 

z  z  0  s  0  0 

It  is  therefore  sufficient  to  consider  only  the  three  integral  equations  (4.  43), 
(4.44)  and  (4.45).  Apart  from  some  slight  simplifications,  these  are  identical  to 
the  equations  obtained  by  Laxpati  (1973)  using  a  scalar  formulation,  and  are  now 
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being  programmed  for  numerical  solution..  Of  the  three  equations,  (4.  43)  and 
(4.46)  are  similar  in  form  and  rather  easy  to  treat.  Indeed,  they  are  not  sig¬ 
nificantly  different  from  those  which  are  employed  in  programs  REST  and  RAM  ID 
(see  Chapter  V).  Unfortunately,  eq.  (4.44)  is  rather  difficult  to  handle  because 
of  the  "second  derivative  singularity"  which  it  contains,  and  an  examination  of 
this  type  of  singularity  is  given  in  Appendix  A. 


V 


COMPUTED  AND  EXPERIMENTAL  RESULTS 

In  the  process  of  assembling  material  for  this  Interim  Report,  we  found 
that  a  great  deal  of  information  had  been  collected  since  the  publication  of  the 
previous  report  (Knott,  Liepa  and  Senior,  1973).  The  theoretical  aspects  of  a 
detailed  examination  of  the  integral  equation  formulation  have  already  been 
presented  in  Chapters  II  through  T  'at  some  of  the  information  we  had  acquired 
seemed  to  fit  in  neither  category,  .v'e  therefore  set  aside  Chapter  V  as  a  reposi¬ 
tory  for  otherwise  unclassifiable  topics. 

Before  presenting  the  details  of  these  assorted  smaller  studies,  we 
first  review  the  salient  events  of  the  predecessor  Contract,  and  of  the  present 
Contract  until  now,  in  Section  5.1.  The  intention  is  to  provide  a  commentary 
that  bridges  the  transition  between  the  two  and  to  show  why  some  of  the  computer 
programs  had  to  be  modified.  The  review  will  refresh  the  memory  of  the  reader 
who  has  seen  the  last  report,  but  will  also  acquaint  the  new  reader  with  previous 
accomplishments.  The  capability  of  the  impedance  boundary  condition  computer 
program  is  assessed  in  Section  5.2  via  comparisons  with  measured  results,  and 
the  effect  of  treating  both  the  leading  and  trailing  edges  of  an  ogival  cylinder  are 
discussed  in  Section  5.  3.  T^e  effect  of  sampling  rate  is  taken  up  in  Section  5. 4, 
along  with  the  demonstration  that,  at  least  for  edge-on  incidence,  it  is  the  lead¬ 
ing  edge  of  the  obstacle  that  dominates  the  scattering  for  E-polarized  incidence. 

5. 1  Summary  of  Previous  Efforts 

Early  in  the  predecessor  Contract  we  were  furnished  a  computer  program 
(RAMI A)  by  AFAL,  based  upon  the  linearized  solutions  of  integral  equations  (2.9) 
and  (2. 17)  derived  in  Chapter  II.  RAM1A  offered  but  two  options:  the  solution  for 
a  circular  cylinder  satisfying  a  (constant)  surface  impedance  boundary  condition 
or  the  solution  for  a  cylinder  of  arbitrary  shape,  but  whose  profile  "ad  to  be 
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specified  point-by-point  as  input  data.  Since  the  investigation  was  to  be  directed 
toward  non-specular  scattering,  circular  cylinders  were  ox  little  interest  and 
the  feeding  of  more  general  profiles  into  the  machine  point-by-point  was  clearly 
an  undesirable  task.  Moreover,  previous  experience  had  shown  that  the  impedance 
must  be  variable  along  the  surface  in  order  to  permit  some  measure  of  control 
over  the  far  field  scattering. 

Since  our  task  required  more  flexibility  than  was  offered  by  RAM1A,  we 
modified  it  so  as  to  more  fully  satisfy  the  existing  needs.  The  profile-generating 
subroutine  was  expanded  and  it  now  constructs  the  body  perimeter  from  a  collection 
of  abutting  straight  line  or  circular  arc  segments  and  assigns  variable  surface 
impedances  along  the  profile  according  to  two  general  mathematical  descriptions. 

A  subroutine  was  also  added  to  the  program  that  prints  out  "quick-look"  plots  of 
the  induced  surface  currents.  To  distinguish  the  program  from  the  original  it 
was  called  RAM  IB. 

RAM1B  was  immediately  tested  for  both  E-  and  H-polarizations  using  a 
conducting  circular  cylinder  as  a  test  obstacle.  The  results  of  these  initial 
test  runs  duplicated  those  of  the  original  program,  and  they  also  agreed  quite 
well  with  the  exact  solution  obtained  by  the  method  of  separation  of  variables 
in  circular  cylindrical  coordinates.  Satisfied  with  the  capability  of  RAM1B, 
we  subsequently  used  it  intensively  to  study  methods  of  reducing  traveling  wave 
returns  from  the  trailing  edge  of  an  ogival  cylinder  and  the  creeping  wave  con¬ 
tributions  from  a  wedge-cylinder  for  H-polarization.  The  results  of  the  study 
showed  that  such  sources  of  non-specular  scattering  can  indeed  be  reduced, 
but  only  at  the  expense  of  treating  a  portion  of  the  surface  which  is  a  significant 
fraction  of  the  wavelength  (Knott  and  Senior,  1973;  Knott,  Liepa  and  Senior, 

1973).  Having  established  favorable  surface  treatments  for  this  polarization, 
we  turned  to  E -polarization,  for  which  the  leading  edges  of  the  bodies  become 
the  dominant  source. 

Preliminary  runs  of  RAM1B  were  for  the  bare  ogival  cylinder  so  as  to 
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establish  reference  levels  against  which  to  gauge  the  success  of  candidate 
surface  treatments.  The  first  test  run  for  the  bare  body  for  E -polarization 
was  unacceptable;  the  surface  fields  rose  to  abnormally  high  levels  at  the 
trailing  edge  and  the  scattered  fields  were  some  10  dB  greater  than  those 
observed  in  experimental  work.  That  the  program  failed  was  corroborated 
by  the  application  of  GTD  as  well,  which,  at  least  within  60  degrees  of  edge-on 
aspects,  agreed  with  the  experimental  patterns  to  within  2  dB.  Clearly,  some¬ 
thing  was  wrong  with  RAMlB  for  E-polarization. 

The  difficulty  had  not  been  resolved  by  the  time  the  Final  Report  of  the 
predecessor  Contract  had  been  prepared,  but  a  series  of  tests  had  been  under¬ 
taken  in  an  effort  to  determine  if  the  problem  could  be  circumvented.  Reasoning 
that  since  the  program  had  worked  quite  well  for  the  circular  cylinder,  it  must 
be  the  sharp  edges  of  the  ogival  cylinder  that  gave  rise  to  the  problem,  so  we  tried 
replacing  the  edges  with  rounded  caps  of  small  radii.  Although  this  ploy  tended 
to  destroy  the  very  feature  that  gives  rise  to  non-specular  scattering,  it  was  at 
once  promising:  the  disagreement  between  measured  and  computed  results 
dropped  from  10  dB  to  only  4  dB.  Then,  reasoning  that  since  the  small  radii 
of  the  caps  actually  resulted  in  dense  surface  sampling  near  the  ends  of  the 
cylinder,  we  verted  back  to  sharp  edges  but  we  clustered  the  sampling  points 
quite  heavily  there.  Subsequent  tests  showed  that  progressively  better  results 
were  obtainable  with  progressively  heavier  sampling  but  this  can  be  carried  only 
so  far;  the  dense  sampling  at  the  edges  greatly  reduces  the  number  of  points 
available  for  distribution  over  the  remainder  of  the  body.  Although  we  were 
able  to  reduce  the  disagreement  between  measured  and  computed  patterns  to 
2  dB,  we  were  disappointed  by  the  devious  means  required  to  do  so. 

In  the  meantime  we  had  received  another  program  (RAMC)  from  AFAL 
which  solved  the  integral  equations  (4. 47)  and  (4. 48)  for  a  metallic  body  in  the 
presence  of  purely  resistive  sheets  for  E -polarization.  Like  RAM1A,  RAMC 
treated  a  specific  geometry,  namely  a  metallic  circular  cylinder  surrounded 
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by  a  concentric  resistive  sheet  of  constant  resistivity,  and  it  was  too  rigid  to  be 
used  in  this  form.  The  program  was  modified  so  as  to  more  fully  suit  our  needs 
and  the  modification  was  named  REST.  Unlike  RAMC,  REST  permits  the  user 
to  place  resistive  sheets  anywhere  (even  within  the  metallic  core! )  and  our  first 
test  of  REST  was  for  a  metallic  ogival  cylinder  with  no  resistive  sheets  at  all. 
When  compared  with  experimental  patterns,  the  computed  data  agreed  quite 
well,  but  it  raised  a  question:  why  did  REST  perform  so  well  while  RAM1B 
did  not?  The  answer  is,  of  course,  that  REST  is  based  on  equations  (4. 47) 
and  (4.48)  for  E -polarization  while  RAM1B  is  based  on  equation  (2. 17),  which 
has  the  second  order  singularity  in  one  term. 

For  the  special  case  of  incidence  normal  to  the  cylinder  axis,  equations 
(2. 16)  or  (2. 17)  may  be  used  for  E-polarization,  of  course,  but  for  numerical 
convenience  (2. 16)  obviously  has  an  advantage.  In  addition,  the  computer 
program  RAM  IB  could  be  greatly  simplified  since  the  generation  of  the  matrix 
elements  for  both  E-  and  H-polarizations  involve  identical  terms.  Thus  we 
put  aside  program  RAMlB  and  constructed  a  far  more  efficient  and  accurate 
version.  It  is  called  RAMD  and  a  listing  of  it  is  given  in  Appendix  B.  As  will 
be  shown  in  a  moment,  RAMD  produces  satisfactory  results  for  both  polarizations, 
at  least  within  the  constraints  of  the  impedance  boundary  condition. 

In  summary,  we  started  the  predecessor  Contract  with  an  impedance 
boundary  condition  computer  program  which  was  then  modified  for  specific 
use  in  non-specular  camouflage  studies  for  H-polarization.  The  present  Contract 
was  started  with  this  program  on  hand,  plus  another  which  solved  the  E-polarized 
resistive  sheet  geometry.  This  second  program  was  also  modified  for  specific 
use  in  E-polarization  studies.  Finally,  a  streamlined  and  compressed  version 
of  the  first  program  was  constructed.  We  now  have  available  programs  REST, 
for  the  resistive  sheet  problem,  and  RAMD  for  the  impedance  boundary  condition. 
As  will  be  pointed  out  presently,  we  are  in  the  process  of  constructing  yet  a  third 
program  that  combines  the  desirable  features  of  both  programs,  valid  for  either 
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polarization.  This  program  will  become,  we  believe,  a  vital  tool  in  the 
remainder  of  our  investigation  of  non-specular  camouflage  treatments. 

5.2  Comparisons  of  Computed  and  Measured  Results 

Programs  REST  and  RAMD  solve  a  two-dimensional  problem,  of  course, 
whereas  laboratory  measurements  can  only  be  made  of  finite,  three-dimensional 
obstacles.  Thus,  in  order  to  compare  measured  and  predicted  patterns,  one  or 
the  other  must  be  scalea  by  the  factor  derived  in  a  previous  report  (Knott,  Liepa 
and  Senior,  1973).  The  relation  between  two-  and  three— dimension  cross 
sections  is 


where  L  is  the  length  of  the  finite  cylinder.  Measurements  were  made  both 
by  AFAL  and  the  Radiation  Laboratory  and,  since  different  frequencies  and 
calibration  units  were  used  at  the  two  installations,  it  is  convenient  to  scale 
the  measured  patterns  down  to  the  two-dimensional  system  instead  of  scaling 
the  computed  patterns  up  to  the  measurements. 

The  AFAL  measurements  were  made  at  7.  5  GHz  for  18— inch  long  cylinders; 
since  AFAL  routinely  presents  its  data  in  dB  relative  to  a  square  meter,  3.  79 
dB  must  be  added  to  the  patterns  in  order  to  obtain  the  effective  two-dimensional 
results.  The  Radiation  Laboratory  measurements  were  made  at  3.  0  GHz  for  a 
cylinder  35  inches  long;  since  our  data  are  routinely  presented  in  dB  relative 
to  a  square  wavelength,  22.  00  dB  must  be  subtracted  in  order  to  carry  out  the 
comparison  with  computed  results.  The  Radiation  Laboratory  patterns  were 
presented  in  the  previous  report,  but  only  one  comparison  was  made  there  with 
predicted  patterns.  The  AFAL  patterns  have  not  yet  been  published. 

Figures  5-1  and  5-2  are  comparisons  of  measured  Radiation  Laboratory 
backscattering  patterns  of  a  metallic  ogival  cylinder  3  wavelengths  wide  for  E- 
and  H-polarizations,  respectively.  Only  90-degree  segments  of  the  patterns 
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are  presented,  and  the  angular  scale  runs  from  edge-on  at  the  left  (zero  degrees) 
to  broadside  incidence  at  the  right  (90  degrees).  Except  for  a  1  dB  error  in  the 
edge-on  region  and  a  2  dB  error  at  broadside,  the  agreement  between  measured 
and  predicted  patterns  for  E-polarization  (.Fig.  5-1)  is  quite  good.  And,  except 
for  the  4  to  5  dB  discrepancy  in  predicting  the  null  near  35  degrees  for  H-polariza- 
tion  (Fig.  5-2),  the  same  kind  of  agreement  prevails.  It  should  be  noted,  however, 
that  the  computed  H -polarized  pattern  seems  to  be  shifted  2  degrees  toward  edge- 
on  incidence.  For  both  polarizations  the  measured  broadside  return  is  less  than 

the  computed  value,  possibly  because  the  target  was  measured  at  a  range  only 

2  • 

half  of  the  customary  2L  /\. 

When  the  ogival  cylinder  is  covered  with  a  layer  of  Emerson  and  Cumming 

SFT-2.  5  absorbing  material  the  agreement  is  not  as  good,  as  may  be  seen  in 

Figures  5-3  arid  5-4.  Based  upon  the  measured  electromagnetic  properties  of 

several  samples  of  this  material  by  AFAL,  a  constant  normalized  surface 

impedance  Z  =  0,  840  +  i0. 187  was  specified  in  the  computer  program  input 
s 

data.  This  impedance  level  corresponds.to  a  normal  incidence  power  reflection 
coefficient  of  -17. 5  dB,  and  if  the  computed  broadside  cross  sections  in  Figures 
5*1  and  5-3  are  compared,  it  will  be  seen  that  a  reduction  of  18.  7  dB  was  pre¬ 
dicted,  reasonably  close  to  that  obtained  from  the  measured  properties  of  the 
material.  However,  a  comparison  of  the  measured  and  computed  patterns  of 
the  coated  ogival  cylinder  for  E-  polarization  show  no  such  agreement.  In  the 
edge-on  region  the  computed  cross  sections  are  of  the  order  of  6  dB  greater 
than  the  measured  values,  and  at  broadside  are  some  5. 5  dB  lower.  Moreover, 
such  pattern  characteristics  as  lcbe  structure  and  null  positions  are  not  well 
produced.  The  disagreement  is  not  quite  so  bad  for  H-polarization  (Fig.  5-4), 
but  the  measured  pattern  does  not  exhibit  the  deep  nulls  predicted  by  the  computer 
program  in  the  broadside  aspect  angle  region.  Interestingly  enough,  the  computed 
broadside  cross  section  reduction  is  only  16. 3  dB  (via  comparison  of  Figs.  5-2 
and  5-4),  and  falls  short  of  the  normal  incidence  value  by  precisely  the  amount 
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that  the  computed  E -polarized  reduction  exceeds  it.  Nonetheless,  we  must 
conclude  that  even  in  the  presence  of  possible  experimental  errors,  the  back- 
scattering  from  the  absorber  covered  ogival  cylinder  is  not  well  modeled  by  a 
constant  surface  impedance  determined  by  the  normal  incidence  properties  of 
the  coating. 

Figures  5-5  and  5-6,  like  Figs.  5-1  and  5-2,  are  for  a  metallic  ogival 
cylinder,  but  in  this  case  the  measured  patterns  were  provided  by  AFAL.  The 
reader  will  note  a  change  in  the  amplitude  scale  on  these  and  the  remaining 
figures,  since  AFAL  uses  a  60  dB  dynamic  range  while  the  Radiation  Laboratory 
uses  40  dB.  In  both  cases  the  cylinder  was  3\  wide  and,  although  the  AFAL 
cylinder  was  shorter  than  that  measured  by  the  Radiation  Laboratory,  the 
difference  in  length  is  automatically  accounted  for  by  the  scaling  factor  men¬ 
tioned  earlier.  The  measured  pattern  for  E -polarization  (Fig.  5-5)  lies  below 
the  computed  pattern  by  about  1  dB  and  were  it  not  for  this  constant  '  Trias",  the 
fit  between  the  two  would  be  remarkable.  The  disagreement  noted  in  the 
Radiation  Laboratory  data  for  H-polarization  (Fig.  5-2)  carries  over  to  the 
case  of  the  AFAL  data  in  Fig.  5-6,  where  again  it  can  be  seen  that  the  computed 
pattern  seems  to  be  shifted  slightly  toward  the  edge-on  aspect  angle. 

Figures  5-7  and  5-8  are  for  the  coated  ogival  cylinder,  but  this  time  the 
coating  was  Emerson  and  Cumming  SFT-11.  0  material.  Based  on  measurements 
of  the  electromagnetic  properties  of  the  layer,  the  normal  incidence  surface  im¬ 
pedance  is  computed  to  be  Z  =0.471  -i0.400,  and  this  was  the  impedance  used 

s 

as  input  data  for  program  RAMD.  Theoretically,  at  least  for  normal  incidence, 
the  radar  cross  section  should  be  reduced  by  6.  97  dB  using  this  impedance,  and 
for  E -polarization  the  computed  reduction  was  in  fact  6,  9  dB,  as  judged  from  a 
comparison  of  Figs.  5-5  and  5-7,  and  for  H-polarization  7.4  dB,  from  Figs. 

5-6  and  5-8.  The  measured  E-polarized  pattern  (Fig.  5-7)  is  about  1.  5  dB 
greater  than  the  computed  one  a(  broadside,  and  from  1.  5  to  2.  0  dB  lower  in  the 
edge-on  region,  but  the  agreement  is  otherwise  quite  good..  For  H-polarization 


fg  2 iiiUM  V  *M»&wUi  ftf  *  H>  '.w 


RELATIVE  POWER  ONE  WAY  (db) 


RELATIVE  FOWER  ONE  WAY  (db) 

-I - - - I _ I _ 1.1:1 _ _ _ I _ I - 1 _ I _ _ _ L L 


(Fig.  5-8),  the  measured  and  computed  patterns  nearly  coincide  at  broadside 
incidence,  but  there  is  some  disagreement  in  the  details  of  the  lobe  structure 
at  edge-on  and  in  the  intermediate  aspect  angle  region. 

Figures  5-9  and  5-10  are  patterns  of  a  metallic  wedge  cylinder  formed 
by  mating  a  25-degree  (total)  angle  wedge  with  a  circular  cylinder  of  electrical 
radius  ka  *  3. 0.  The  measured  pattern  for  E -polarization  (Fig.  5-9),  like  that 
of  the  metallic  ogival  cylinder,  lies  about  1.  5  dJB  below  the  computed  one  and 
were  it  not  for  this  fairly  constant  difference,  the  two  patterns  would  be  in 
almost  perfect  agreement.  The  agreement  for  H-polarization  (Fig.  5-10)  is 
not  quite  as  good,  with  some  of  the  computed  nulls  shifted  toward  the  edge-on 
incidence  side  of  the  pattern  and,  in  particular,  with  a  6  dB  disagreement  in 
amplitude  at  the  edge-on  aspect.  Note  that  since  the  wedge  angles  of  both  the 
wedge-cylinder  and  the  ogival  cylinder  are  the  same,  the  edge-on  E-polarized 
returns  of  both  bodies  are  the  same.  Figures  5-1,  5-2,  5-5,  5-6,  5-9  and  5-10 
all  suggest  that  the  computer  program  performs  quite  well  for  metallic  bodies 
for  E -polarization  and  that,  while  less  accurate,  it  may  still  produce  acceptable 
patterns  for  H-polarization.  Some  of  the  differences  noted  may  be  due  to  align¬ 
ment  difficulties  and  near  field  effects,  and  not  least,  to  the  neglect  of  end 
effects  in  scaling  the  three-dimensional  patterns  down  to  the  two-dimensional 
-jomputed  results. 

Figures  5-11  and  5-12  display  the  wedge-cylinder  patterns  when  the 
body  is  coated  with  the  SF-11.0  material  mentioned  earlier.  The  agreement 
for  E-polarization  (Fig.  5-11)  ranges  from  1.  5dB  at  the  speculai  angle  to  3  dB 
or  so  in  the  regions  of  the  nulls,  and  is  as  good  as  was  noted  for  the  ogival 
cylinder.  The  discrepancy  for  H-polarization  (Fig.  5-12)  is  somewhat  greater, 
about  3dB,  at  the  specular  angle,  but  rises  to  as  much  as  10  dB  in  the  edge-on 
region..  Comparisons  of  the  specular  cross  sections  in  these  figures  show 
that  the  computed  reduction  is  6.:  ldB  for  E-polarization,  which  is  0.  9  dB 
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short  of  what  might  be  expected  for  n.'mnal  incidence  as  judged  from  the 
material  properties,  and  is  7. 5  dB  for  H-polarization,  which  is  0. 6  dB 
better.  The  measured  cross  section  reductions  considerably  exceed  the  com¬ 
puted  ones,  however. 

The  agreement  between  measured  and  computed  patterns  is  quite  good 
for  metallic  bodies,  nearly  as  good  for  coated  bodies  for  E-polarizatio  and 
perhaps  disappointing  for  coated  bodies  for  H-polarization.  We  believe  the 
errors  are  traceable,  in  part,  to  the  shortcomings  of  the  surface  impedance 
boundary  condition,  which,  while  conceivably  justifiable  over  smooth  surfaces 
of  large  radii  of  curvature  for  sufficiently  lossy  materials,  is  not  necessarily 
valid  near  edges.  From  a  study  of  the  integral  equations  presented  in  Chapter 
II,  one  would  conclude  from  that  a  surface  impedance  boundary  condition  tends 
to  irterchange  the  roles  of  E-  and  H-polarized  scattering  sources.  Hence  the 
surface  impedance  would,  for  example,  increase  the  trailing  edge  return  for 
E -polarization.  If  this  were  truly  the  case,  the  design  might  have  to  aim  for 
a  surface  impedance  of  about  unity,  since  this  would  be  optimum  for  arbitrary 
polarizations. 

However,  the  experimental  data  in  Figs.  5-1  through  5-12  show  no  such 
inclination.  In  every  case  the  radar  cross  section  in  the  edge-on  region  was 
reduced;  no  enhancement  took  place,  as  might  have  been  predicted  from  the 
duality  of  the  impedance  boundary  condition  integral  equations.  Therefore  the 
effective  surface  impedance  near  the  edges  v'as  different  for  the  two  polarizations, 
even  though  the  edges  were  treated  with  the  same  physical  materials.  This  is 
indeed  fortunate,  for  it  gives  us  the  kind  of  control  we  require:  the  optimum 
surface  impedance  for  one  polarization  is  mani'estly  not  optimum  for  the  other. 
Thus,  instead  of  being  disappointed  in  the  lack  of  agreement  between  the 
measured  and  predicted  patterns  of  these  coated  objects,  we  take  heart:  the 
disagreement  shows  that  selective  control  is  still  available. 
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5.3  Single  and  Double  Edge  Treatment  of  the  Ogival  Cylinder 

Although  it  is  doubtful  that  an  impedance  can  be  defined  near  the  edges 
of  a  coated  body  on  the  basis  of  the  normal  incidence  properties  of  the  coating 
alone,  there  is  no  question  that  some  impedance  can  be  assigned.  What  is  net 
clear  is  how  we  might  deduce  this  impedance  from  a  knowledge  of  the  material 
properties  and  the  local  edge  geometry,  and  indeed,  it  is  via  the  expansion  of 
program  TWOD  that  we  hope  to  answer  the  question  in  the  coming  months. 

The  impedance  boundary  condition  integral  equations  do  not  in  themselves 
provide  the  answer.  They  merely  offer  the  promise  of  providing  a  solution  to  a 
given  scattering  problem  for  whatever  particular  impedance  is  selected  to 
represent  the  relation  between  the  tangential  surface  fields.  Not  knowing  at 
present  how  the  H-polarized  impedance  near  an  edge  is  related  to  the  E- 
polarized  impedance  there,  we  shall  assume  for  the  time  being  that  they 
are  equal.  We  now  pose  the  question,  how  would  an  optimum  treatment  for 
one  polarization  affect  the  performance  for  the  other,  assuming  that  the 


impedance  for  the  two  polarizations  is  the  same  ? 

To  obtain  some  idea  of  the  possible  degradation  in  performance,  an 

ogival  cylinder  was  chosen  as  a  test  obstacle  because  the  offending  source  of 

non-specular  scattering  shafts  from  the  leading  edge  to  the  trailing  edge  as 

the  polarization  is  rolled  from  E  to  H.  We  assume  that  a  parabolic  (square 

law)  impedance  distribution,  rising  to  a  terminal  value  of  Z  =  2.  O  +  iO  at  the 

s 

edge,  is  optimum  for  H-polarization  when  1  \  of  the  trailing  edge  surface  is 
treated,  the  constraints  of  the  optimization  being  due  to,  say,  external  con¬ 
siderations  such  as  weight,  layer  thickness,  and  the  like.  Normally,  only 


the  trailing  edge  should  be  treated,  but  if  the  cylinder  must  also  be  viewed 
fi’om  the  rear,  then  the  leading  edge  should  be  included. 

Figure  5-13  summarizes  the  results.  As  seen  on  the  left  side  of  the 
diagram,  a  trailing  edge  treatment  alone  (the  alternating  solid-dashed  line) 
reduces  the  return  by  2  dB  at  precisely  edge-on  incidence,  but  out  to  about 
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COMPARISON  OF  SINGLE  AND  DOUBLE  EDGE  TREATMENTS  FOR  H-POLARIZATION 
USING  SURFACE  IMPEDANCES. 


40  degrees  an  average  of  10  dB  or  so  is  available.  The  treated  edge  becomes 
a  "leading"  edge  if  we  shift  from  the  left  to  the  right  side  of  the  diagram,  and 
there  it  can  be  seen  that  the  cross  section  is  actually  enhanced  oy  5  dB  over 
the  bare  body  level,  that  a  deep  null  is  created  at  165  degrees  (corresponding 
to  15  degrees  from  edge-on),  and  that  the  averaged  cross  section  values  in  the 
140  to  180  degree  range  have  not  been  materially  affected.  In  fact,  an  envelope 
sketched  through  the  peaks  of  the  lobes  would  seem  to  describe  the  peak  ampli¬ 
tude  of  the  bare  body  echo  as  well  as  that  of  the  body  with  a  single  treated  edge 
(when  that  edge  is  presented  toward  the  observer).  Thus  a  favorable  trailing 
edge  treatment  is  completely  sacrificed  if  the  edge  is  viewed  from  the  rear. 

This  being  so,  it  might  be  expected  that  treating  both  edges  of  the  body 
would  be  useless,  for  then  a  coated  edge  would  be  presented  to  the  observer 
regardless  of  whether  he  viewed  the  cylinder  from  the  front  or  the  rear.  Such 
is  not  the  case,  however,  as  can  be  seen  from  the  level  of  the  dashed  line  in 
Fig.  5-13,  representing  the  effect  of  a  double  edged  treatment.  In  the  inter¬ 
mediate  aspect  angle  range  from  45  to  135  degrees,  the  double  edged  treatment 
smoothes  out  the  peaks  and  nulls  of  the  bare  body  pattern,  and  the  averaged 
cross  section  levels  are  scarcely  different  from  those  of  the  metallic  objects. 
But  in  the  edge-on  regions,  a  7  dB  reduction  is  still  available,  although  there 
is  a  substantial  enhancement  of  about  6  dB  at  precisely  edge-on.  We  conclude 
that,  although  a  trailing  edge  treatment  offers  little  benefit  if  viewed  from  the 
wrong  direction,  a  double  edged  treatment  helps  to  recover  the  loss  in  perfor¬ 
mance.  The  recovery  is  not  complete,  however,  and  some  performance  is 
sacrificed. 

This  does  not  seem  to  carry  over  to  the  case  of  E-polarization,  primarily 
because  the  trailing  edge  does  not  appear  to  be  a  source  of  scattering.  As 
shown  in  Fig.  5-14,  for  example,  nothing  is  lost  by  treating  both  edges  of 
the  cylinder.  Aside  from  a  shift  in  the  positions  of  the  peaks  and  nulls  of  the 
pattern,  the  averaged  amplitudes  of  the  scattering  from  both  the  single  and 
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USING  SURFACE  IMPEDANCES. 


double  edged  treatments  are  about  the  same,  but  the  double  edged  version  does 
have  the  virtue  of  reducing  the  return  in  both  the  forward  and  backward  directions. 
The  cross  section  reduction  available  for  E-polarization  is  substantial,  amounting 
to  nearly  20  dB  within  a  40-degree  range  of  edge-on  incidence.  Thus,  at  least 
under  the  "optimized"  conditions  specified  for  this  particular  case,  H-polarization 
is  likely  to  suffer  most  from  the  double  edged  treatment,  with  the  E-polarized 
cross  section  being  scarcely  of  concern. 

The  above  results  were  obtained  using  program  RAMD  and  it  was  of 
interest  to  perform  the  same  kind  of  test  using  program  REST,  in  which 
resistive  sheets  are  employed.  We  chose  a  sheet  1\  wide  having  a  parabolic 
resistance  variation,  with  the  terminal  (outer  edge)  resistance  equal  to  twice 
the  impedance  of  free  space.  Thus  the  sheet  width  and  the  resistance  variation 
were  analogous  to  the  surface  impedance  treatment  used  in  conjunction  with 
RAMD.  The  results  are  displayed  in  Fig.  5-15  and,  like  those  in  5-14,  the 
double  edged  treatment  entails  little  or  no  loss  in  performance  compared  with 
that  obtained  by  treating  a  single  edge. 

There  are  two  minor  differences  between  the  results  of  the  resistive 
sheet  and  surface  loading  methods,  however.  The  first  is  that  the  resistive 
sheet,  since  it  projects  outward  from  the  bare  body,  increases  the  effective 
length  of  the  scatterer.  Consequently  the  lobe  structure  tends  to  be  more 
detailed  and  the  broadside  cross  section  is  slightly  greater.  For  the  1\  sheet 
width  used  in  the  test  of  Fig.  5-15,  this  enhancement  amounts  to  about  1.5  dB. 

The  second  difference  is  that  the  impedance  loading  method,  when  extended 
over  a  surface  distance  equal  to  the  sheet  width  in  the  above  test,  offers  a 
slightly  better  cross  section  reduction.  The  improvement  amounts  to  2  to  4 
dB  in  the  end-on  region  and  it  persists  into  the  broadside  region  where  the  resis¬ 
tive  sheet  actually  produces  an  enhancement. 

As  was  done  under  the  predecessor  Contract  for  the  trailing  edge  with 
H-polarization,  we  also  studied  the  effect  of  the  impedance  distribution  near  the 
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leading  edge  for  E-polarization  in  the  present  Contract.  The  study  was  not 
exhaustive,  but  the  results,  some  of  which  are  shown  in  Fig.  5-16,  suggest 
that  the  leading  edge  treatment  is  not  nearly  as  difficult  to  perform  for  E- 
polarization  as  was  the  ti*ailirg  edge  loading  for  H-polarization.  The  angle 
of  incidence  was  held  fixed  at  25  degrees  from  edge-on  (as  in  the  H -polarized 
studies  previously  mentioned)  and  two  types  of  distributions  were  imposed. 

The  length  of  the  treated  surface  ranged  from  0.  5X  to  1. 5X,  the  latter  covering 
half  the  body. 

The  four  curves  plotted  in  Fig.  5-16  are  typical  of  all  the  data  collected 

for  this  case.  As  the  figure  shows,  slightly  less  performance  is  obtained  for 

shorter  surface  loadings,  and  it  seems  to  matter  little  whether  a  linear  or  a 

parabolic  distribution  is  employed.  The  maximum  impedance  used  in  the 

study  was  Z  =  2.  0  +  i0  (at  the  leading  edge)  but  the  curves  seem  to  be  continuing 
s 

downward  monotonically.  They  suggest  that  the  scattering  does  indeed  arise 
from  the  edge  itself  and  that  the  loading  might  be  compressible  into  a  narrow 
strip  of  surface  near  the  edge  if  the  impedance  can  be  taken  to  a  high  enough 
level.  A  cross  section  reduction  of  13  dB  is  apparently  attainable  if  the  leading 
edge  impedance  is  as  low  as  1.  3  times  that  of  free  space.  The  form  of  the 
curves  is  much  like  that  provided  by  resistive  sheets  placed  in  front  of  the 
cylinder. 

5.4  Reducing  the  Number  of  Sampling  Points 

The  expansion  of  program  TWOD  to  include  magnetic  was  well  as  dielectric 
properties  of  physical  materials,  for  H-polarization  as  well  as  E,  will  produce 
a  useful  tool  for  our  task.  The  program  will  be  based  upon  the  integral  equations 
for  thin  layers,  but  it  is  not  inconceivable  that  relatively  thick  coatings  may 
have  to  be  studied,  or  even  coatings  whose  thicknesses  are  variable.  Hence 
it  may  be  necessary  to  synthesize  such  coatings  by  stacking  up  several  thin 
layers.  This  amounts  to  sampling  in  volume  instead  of  area  and,  unless  the 
total  number  of  sampling  points  is  increased,  it  will  require  that  the  points  be 
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redistributed.  Some  parts  of  the  body  will  therefore  be  sampled  le  -s  densely 
than  others,  and  we  ran  a  sequence  of  tests  to  establish  how  the  ccmouted  far 
fields  might  be  affected  by  sparse  sampling. 

Again  the  ogival  cylinder  was  chosen  as  a  test  obstacle  and  program 
REST  was  used  to  compute  the  backscattering.  No  resistive  sheets  were 
deployed,  hence  the  patterns  produced  were  of  the  bare  body.  The  total  number 
of  sampling  points  on  the  profile  was  increased  in  steps  of  six  from  a  minimum 
of  24  to  a  maximum  of  100,  producing  nominal  surface  sampling  rates  from  4 
to  16%  per  wavelength.  Assuming  that  the  highest  sampling  rate  produces 
"perfect"  results,  the  errors  produced  by  less  frequent  sampling  can  be  estimated 
by  comparing  the  computed  backscattering  patterns.  This  comparison  is  shown 
in  Fig.  5— IV  for  a  selection  of  5  sampling  rates. 

The  errors  are  greatest  in  the  intermediate  aspect  angle  range  and  are 
smaller  at  both  the  edge-on  and  broadside  aspects.  Expectedly,  ‘he  error 
increases  with  decreasing  sampling  rates.  From  the  data  display »  one  finds 
that  10  samples  per  wavelength  yield  0.  5  dB  accuracy  for  the  entire  aspect 
angle  range,  but  if  only  edge-on  incidence  is  of  interest,  as  few  as  samples 
per  wavelength  will  suffice.  Figure  5-18  is  a  plot  of  the  errors  incurred  at 
three  discrete  aspect  angles,  one  of  the  angles  (0  =  58  degrees/  being  that  at 
which  the  error  is  greatest.  It  should  be  emphasized  that  these  p.'ots  are  for 
a  bare  ogival  cylinder  and  that  the  behavior  depicted  may  not  necessarily  apply 
to  other  bodies  or,  indeed,  to  the  same  body  when  coated. 

If  it  is  only  the  leading  edge  that  dominates  the  far  field  scattering,  as  is 
true  for  E -polarization,  it  may  be  possible  to  use  a  small  wedge-cylinder,  as 
suggested  in  Fig.  5-19.  In  this  case,  the  sampling  rate  was  held  fixed  at  about 
16  samples  per  wavelength,  and  the  electrical  size  of  the  body  was  varied  from 
ka  =  0. 625  to  ka  =  3.  0.  The  number  of  sampling  points  ranged  from  20  to  96, 
yet  the  edge -on  return  varied  by  less  than  0. 3  dB  throughout  the  range.  These 
data  suggest  that  a  great  number  of  sampling  points  may  be  released  for  other 
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FIG.  5-17:  ERRORS  IN  COMPUTED  FAR  FIELD  PATTERN  OF  AN 
OGIVAL  CYLINDER  AS  A  FUNCTION  OF  ASPECT 
ANGLE. 
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FIG.  5-18:  ERRORS  IN  COMPUTED  BACKSCATTERING  PATTERN  OF  AN 
OGIVAL  CYLINDER  AS  A  FUNCTION  OF  SAMPLING  RATE. 
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FIG.  5-19:.  EDGE -ON  RETURN  OF  A  METALLIC  WEDGE -CYLINDER 
FOR  E-POLARIZATION. 
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parts  of  the  body  (e.  g.,  for  the  coating),  provided  the  original  scattering 
source  is  preserved.  Since  the  dominant  source  then  is  the  leading  edge, 
the  size  of  the  cylinder  is  relatively  unimportant,  but  in  other  cases,  such 
as  H-polarization  for  the  same  body,  this  is  no  longer  true.  Thus,  in  the 
event  that  it  becc  .es  necessary  to  re-distribute  sampling  points,  there  are 
certain  instances  for  which  it  can  be  accomplished  with  but  a  small  sacrifice 
in  accuracy.  The  re-allocation  is  not  possible  in  all  cases,  of  course,  and  we 
may  have  to  search  for  other  techniques  should  the  need  arise. 


75 


VI 


CONCLUSIONS 

In  this  report  we  have  described  most  of  the  work  carried  out  since  the 
publication  of  the  Final  Report  (Knott  et  al.,  1973)  uuuer  the  predecessor  Con¬ 
tract.  One  of  the  key  features  of  this  study  is  the  detailed  examination  of  the 
integral  equations  in  two-dimensional  scattering  presented  in  Chapters  n  and 
IV.  For  an  impedance  boundary  condition  imposed  at  the  surface,  the  analysis 
has  revealed  that  there  is  a  choice  of  integral  equations  appropriate  to  a  solid 
body,  and  by  using  this  fact  we  were  able  to  circumvent  the  difficulties  previously 
encountered  using  program  RAMIA(B)  for  E-polarization.  The  result,  now 
designated  RAMD,  was  a  far  more  efficient  and  effective  program  and  is  listed 
in  Appendix  B.  RAMD  has  been  used  to  examine  the  effect  of  tapered  impedances 
on  the  scattering  from  ogival  and  wedge-cylinders  for  E-polarization.  Some  of 
the  data  obtained  are  presented  in  Chapter  V,  and  are  compared  with  experi¬ 
mental  data. 

One  of  the  shortcomings  of  an  impedance  boundary  condition  is  the  difficulty 
of  relating  the  surface  impedance  to  the  electromagnetic  properties  of  a  material 
necessary  to  simulate  it.  To  provide  a  more  explicit  connection  to  material 
properties,  we  have  given  considerable  attention  to  the  concept  of  resistive  sheets 
used  either  individually,  or  in  combination  to  simulate  a  layer  of  finite  thickness. 
For  the  problem  of  (electrically)  resistive  sheets  with  E-polarization,  a  computer 
program  is  available  and  is  designated  REST.  We  have  employed  it  extensively 
and  successfully  to  assess  the  sheet  performance,  but  the  results  obtained  are 
reserved  for  a  future  report.  However,  it  is  important  to  examine  the  behavior 
of  these  same  sheets  for  H-polarization,  and  since  we  can  also  conceive  of 
'magneticallv  resistive'  sheets,  to  include  these  in  our  study  as  well. 

The  derivation  of  the  integral  equations  for  both  types  of  resistive  sheets, 
either  in  isolation  or  in  the  presence  of  a  body  having  an  impedance  boundary 
condition  imposed  at  its  surface,  is  presented  in  Chapter  IV.  In  essence,  this 
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serves  to  set  the  stage  for  the  work  we  are  now  doing,  and  the  construction  of 
a  computer  program  for  generalized  resistive  sheets  is  now  almost  complete. 
Although  at  least  one  of  the  integral  equations  is  characterized  by  a  'second 
derivative  singularity'  which  is  difficult  to  handle,  the  analysis  and  testing 
described  in  Appendix  A  shows  how  it  can  be  treated  numerically.  The 
fruits  of  this  work  will  be  the  subject  of  the  next  report. 
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APPENDIX  A 

SOME  COMPUTATIONAL  CONSIDERATIONS 


Many  of  the  integral  equations  derived  in  Chapters  II  and  IV  resulted  from 
a  limiting  operation  of  the  form 


lim  ,__3_ 
£->  C  3n 


f(s')HQ^(kr)ds' 

^C 


(A.  1) 


if  £  is  not  on  C,  the  differentiation  and  integration  can  be  interchanged  to  give 


where 


lim 

£-^C 


f(s*)K(s,  s')ds' 

•c 


K(s, s*)  =  k^-rlH^V)  -  -ktf-rjH^fer)  , 


(A.  2) 


(A.  3) 


but  in  the  limit  as  the  observation  point  approaches  C,  the  kernel  is  characterized 
by  a  "first  derivative  (non-integrable)  singularity"  and  is  infinite  at  s’  =  s.  How¬ 
ever,  (A.  2)  and,  hence,  (A.  1)  can  be  evaluated  by  treating  analytically  a  neigh¬ 
borhood  of  the  point  s'  =  s,  and  when  this  is  done  we  find 


lini 
£  — ^  C 


f(s*)(n-r)H^  (kr)d(ks’)  =  ±2if(s)+  -V  f(s')(n  r)H^  (kr)d(ks') 


(A.  4) 


with  the  upper  or  lower  sign  according  as  £  approaches  C  in  the  direction  t  n, 
respectively.  The  slash  across  the  integral  sign  denotes  the  Cauchy  principal 
value,  and  since  the  first  term  is  simply  the  self  cell  contribution,  the  right  hand 
side  of  (A.  4)  is  quite  convenient  for  computation. 
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Unfortunately,  several  of  the  key  integral  equations  involve  a  second 


derivative  kernel  in  addition  or  instead,  for  example 


lim 

£-»C 


f(s')  sOds* 

on 


'C 


(A,  5) 


where  K(s,s!)  is  given  in  (A.  3) ,  and  this  is  not  quite  so  easy  to  handle.  Even 
the  existence  of  the  limit  is  mathematically  in  question,  but  the  form  that  preceded 
(A.  5) ,  viz. 


lim 

£~>C 


_a 

On 


f(s')K(s,  s’Jds* 

b 


(A.  6) 


is  capable  of  justification  and  physically  represents  the  field  of  a  current  sheet 
evaluated  at  its  surface. 

To  see  how  we  can  treat  it,  consider 


i(x,y) 


(A.  7) 


appropriate  to  a  planar  sheet.  If  x  is  on  C  and  A  is  a  small  interval  of  width 
26  centered  on  x'  =  x,  l£x,y)  can  be  written  as 

I(x,y)  =  I1(x,y)+I2(x,y) 

where 

^(x.y)  =  ^  ^  J(x*)  ~~2  Ho^Ck  \Z^-x’)2+  y^dx*  ,  (A. 8) 

JC-A 

I2(x,y)  =  £  \  J(x')H^  ^k'J(x-x,)2+y2')dx»  .  (A.  9) 

*  Ja 
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Taking  1^  first  and  using  the  differential  equation  for  the  Hankei  function,  we 


Ix0c,y)  = 


Jfc')  ~yT  ^'(k^fe-x'^+y4, 

WCx-x’)  +y 


,  2(1)/  r  Z2~~2\  dx» 

-ky  Hq  (kV(x-x')  +y  ft-  -V- 

J  (x  -x1)  +y 


and  since  the  integral  is  convergent  for  all  y  including  y  =  0, 


urn  ,  .  . 

y-»o  1i(x'y)  -- 


H^Vlx-x'l) 

M  dx'  • 


(A.  10) 


2  .  2 


Turning  now  to  we  remark  that  if  k  V<5  +y  «  1  and  J(xT)  is  slowly  varying 
in  the  vicinity  of  x*  =  x, 


J (x')Hq1  '  (k\J (x  - x’)2  +y2 


dx'  ~  —  J(x)\  logf~  ^/(x-xO^+y2  )dx* 


=  f  J(x)|2i6|log  (Y^+y2)-^  -ylog^“ 


where  T=  eT  with  y  =  0,5772. . .  (Euler's  constant).  On  carrying  out  the  dif¬ 
ferentiation  with  respect  to  y. 


I  (x,y)  -  -  - - - —  J (x) 

ffk(6  +y  ) 


from  which  we  obtain 


,!£o  «**•»»  *-&*»  ■ 


(A.  11) 


and  hence 
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iim  4i  .  . 

I(x,y)  —  -  — —  J(x)  - 
y  — 7rk6 


H^(kjx-x'j) 

J(x')  - n —  dx* 


(A.  12) 


It  car.  b°  verified  that  this  result  is  unaffected  if  one  or  both  of  the  differentiations 
are  carried  out  prior  to  integration. 

The  right  hand  side  of  (A.  12)  is  not  inconvenient  for  numerical  evaluation, 

but  in  contrast  to  (A.  4),  there  is  no  guarantee  of  increased  accuracy  as  the  cell 

size  26  is  reduced.  The  first  term  is  the  self  cell  contribution  and  is  inversely 

proportional  to  the  cell  size.  The  factor  multiplying  J(x*)  in  the  integrand 
_2 

increases  as  (x-x‘)  as  x'  approaches  x,  but  since  the  self  cell  is  excluded,  the 
integrand  is  finite  for  all  x'.  However,  over  the  immediately  adjacent  cells  the 
integrand  can  vary  by  as  much  as  an  order  of  magnitude,  and  it  will  almost  cer¬ 
tainly  be  necessary  to  take  account  of  this  variation  through  subdivision  of  these 
cells  or  by  integration  over  them. 

This  approach  can  be  tested  using  the  simple  example  of  a  plane, 
electrically  resistive  sheet  illuminated  by  an  H-polarized  plane  wave  at  non¬ 
glancing  incidence.  If  the  sheet  occupies  the  portion  0  <  x_<L  of  the  plane  y  -  0 


,ri  a  -ik(xcosa+ysinQ') 
H  =  z  e 


(A.  13) 


eq.  (A.  12)  can  be  written  as 


-ikxcosc?  \n  \  lim  1  9 

smoe  =  YJR(x)J(x)  -  .  _  — - ~ 

0  y  — >  0  4k  g  2 


j(xf)H' 


'(kj/fe- 


,.2  2  V  , 

xf)  +y  Jdx* 


and  hence,  from  (A.  12), 


-ikx  cos  a  _  f  .  .  i  1  T/  .  ,1 

“  [Y0RW+SsJJW+4 


H^l\kjx-x'j) 


dx*  . 


(A.  14) 
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Equation  (A.  14)  has  been  programmed  for  solution  and  the  current  J(x) 
computed  in  the  particular  case  of  a  highly  conducting  (R  =  0.  01  ohm)  strip  3A 
long  for  plane  wave  incidence  at  the  angle  a  =  155  degrees,  corresponding  to  the 
peak  of  the  traveling  wave  lobe  in  the  backscattering  pattern.  From  an  examina¬ 
tion  of  the  data,  it  is  believed  that  48  sampling  points  (so  that  26  =  A/ 16, 
implying  k6  =  0. 196),  with  the  cells  immediately  adjacent  to  the  self  cell  sub¬ 
divided  into  25,  and  the  next  5  cells  on  either  side  each  subdivided  into  10,  are 
sufficient  to  give  data  which  are  effectively  ''exact".  The  resulting  curve  for 
|j(x)|  is  s..  wn  in  Fig.  A-i  and  has  the  character  expected  in  a  traveling  wave 
situation.  The  effect  of  reducing  the  number  of  subdivisions  is  illustrated  in 
Table  A-l,  from  which  it  appears  that  if  we  are  to  keep  the  errors  less  than  1  per¬ 
cent,  it  is  necessary  to  subdivide  two  cells  on  either  side  of  the  self  cell,  with  the 
nearer  being  more  finely  subdivided.  Many  of  these  same  tests  have  also  been 
performed  with  24  and  32  sampling  points  (or  cells).  It  is  found  that  the  errors 


TABLE  2-1;  ERRORS  IN  jj  | 


subdivision 

aver.  (°/o) 

max.  (°/o) 

1  x  25,  5  x  10 

as  sun 

led  0 

1  x  25,  4  x  10 

0.06 

o.os 

1  x  25,  3  x  10 

0. 11 

0.  18 

1  x  25,  2  x  10 

0.09 

0.21 

1  x  25,  l  x  10 

0.32 

0.61 

1  x  10,  1  x  10 

1.53 

2.23 

1  x  25 

3.39 

4,  85 

1  x  10 

i  4.53 

6.38 

are  significantly  increased.  As  an  example,  with  N  =  24  and  six  cells  on  each 
side  subdivided  as  in  the  first  line  of  Table  A-l,  the  average  and  maximum  errors 
are  7.31  and  16. 03  percent  respectively. 

In  any  subdivision  process  the  integrand  must  be  computed  a  corresponding 
number  of  times,  and  any  massive  subdivision  is  therefore  an  inefficient  and  time 
consuming  procedure.  If  the  strip  is  divided  into  48  cells,  subdivision  of  the 
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FIG  A_1;  TOTAL  CURRENT  ON  A  HIGHLY  CONDUCTING  (R  =  0,01  n)  STRIP  3X  WIDE  ILLUMINATED 
BY  AN  H-POLARIZED  PLANE  WAVE  WITH  a  =  155°. 


ceils  adjacent  and  next-to-adjaeent  to  the  self  cell  by  25  and  10  respectively 
requires  almost  1500  additional  computations  of  the  Hankel  function  even  when 
the  symmetry  of  the  matrix  is  taken  into  account,  and  this  increases  consider¬ 
ably  the  time  involved  in  generating  the  matrix  elements.  However,  subdivision 
is  just  a  simple  method  of  numerical  integration,  and  since  the  unknown  current 
J(x)  is  still  given  its  value  at  the  midpoint  of  -he  cell,-  the  function  being  integrated 
is  entirely  known.  This  suggests  the  possibility  of  an  analytical  evaluation.  The 
contribution  of  the  kernel  from  the  nth  ceil  on  either  side  of  the  self  cell  in  eq. 

(A.  14)  .. 

r<2"+1){  Hu)to) 

I  =  \  - - - dt  (A.  15) 

n  1  t 

J(2n-1)6 

and  from  the  recurrence  relations  for  the  Hankel  functions  (or,  alternatively, 
from  the  differential  equation  for  H^(kt)  ),  we  have 

pk6 

1  =  H  jl)  ([2n  -  l]  k6^  -  H^  ([2n  + 1]  kb)  +  \  H^(2nk6+t)  dt  .  (A.  .6) 

J-k6 

The  integrand  that  now  remains  is  much  more  slowly  varying  than  the  one  in 
(A.  15)  ,  and  it  is  not  unreasonable  to  approximate  (A..  16)  by 

1  ~  H!1)  ( [2n  ~  J]  k6)  "  Hii) ([2n  + 1] k6)  +  2k 6 H^ (2nk6)  , 

Somewhat  surprisingly,  when  this  formula  is  used  for  n  =  1,2,  i.e,,  for 
two  cells  on  either  side  of  the  self  cell,  with  k6  =  0.  196,  the  current  is  in  error 
by  more  than  10  percent.  Most  of  this  error  is  attributable  to  the  evaluation  of 
the  integral  in  (A.  16) ,  and  is  not  substantially  reduced  when  the  Hankel  function 
H^(x)  is  replaced  by  its  iogaiitbmic  approximation  for  small  arguments  and 
the  integration  carried  out  analytically.  On  the  other  hand,  when  the  cells  cor¬ 
responding  to  n  =  1,  2  are  subdivided  into  5  and  3  parts  respectively,  the 
average  and  maximum  errors  drop  to  0.:19  and  0.42  percent  respectively,  and 
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even  with  the  coarser  subdivision  3  and  2,  the  errors  are  still  only  0.  37  and 
0.G8  percent.  Such  a  performance  is  comparable  to  that  achieved  on  dividing., 
the  same  cells  in  (A.  15)  into  25  and  10  parts  (see  line  5  of  Table  A-i),  and 
requiies  only  about  300  additional  computations  of  the  Hankel  functions  in  place 
of  1500.  As  an  alternative  to  subdivision,  integration  using  the  three-point 
Simpson  formula  has  been  tried.  When  applied  to  the  cells  corresponding  to  n  =  1 
and  2  in  eq.  (A.  16) ,  it  turns  out  to  be  just  slightly  more  efficient  than  subdivision 
of  the  same  cells  into  3  and  2  parts  respectively,  but  produces  much  better 
accuracy,  the  average  and  maximum  errors  in  |j  |  being  0.  12  and  0  40  percent. 

In  contrast,  when  the  same  integration  formula  was  applied  to  the  original  expres¬ 
sion  (A.  15)  for  the  cell  contributions,  the  errors  were  larger  by  two  orders  of 
magnitude. 

All  of  the  above  constitutes  the  direct  approach  to  an  integral  of  the  form 
shown  in  eq.  (A.  7) ,  and  though  it  has  turned  out  to  be  relatively  straightforward, 
it  is  not  the  only  way  in  which  the  second  derivative  singularity  can  b*>  handled. 

If,  for  example,  we  gc  back  to  (A.  7)  and  use  the  fact  that 


the  normal  derivatives  can  be  replaced  by  tangential  derivatives  and  one  (or  both) 
of  them  eliminated  by  partial  integration.  This  method  has  been  explored  in  some 
detail.  Integrating  only  once  by  parts,  we  obtain  an  integral  expression  for 
I(x,y)  involving  the  current  and  its  first  derivative,  plus  contributions  from  the 
end  points  of  the  range  of  integration.  The  more  singular  integral  is  that  con- 

0j 

taining  -^7  ,  but  in  the  limit  as  y-*0  the  self  cell  contribution  can  be  evaluated 
analytically  and  turns  out  to  be  zero.  Hence,  for  0  x  <-  L, 
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II 


y->0 


iOc.y)  = 


\  /  1  \ 

-k  \  j(x’)H  7(k  jx-x'J  )dx' 


+J(L)H,  (k  [L-x])+J(Q>H.  (kx) 
1  W  "  1 


(A.;  17) 


and  though  we  are  now  left  with  twice  as  many  unknowns  as  we  had  before,  the 

assumption  of  a  locally  linear  variation  of  J(x*)  enables  us  to  express  — —  in 

ax’ 

terms  of  the  currents  over  the  neighboring  cells. 

Unfortunately,  the  end-point  contributions  sose  a  problem.  Each  is 
infinite  when  the  observation  point  is  at  the-  end  in  question,  and  while  it  can  be 
argued  that  in  the  moment  method  the  observation  never  does  lie  actually  at  an 
end,  the  contribution  will  still  increase  indesandfcely  with  decreasing  cell  size  In 
effect,  we  now  have  a  self  cell  contribution  analogous  jet  oA.  11),  Qut  appropriate 
only  when  the  self  cell  is  at  an  end.  If,  on  the  osaer  hand,  the  self  cell  contribu¬ 
tion  is  separated  out  prior  to  partial  integration,  the  result  obtained  is  no  dif¬ 
ferent  from  integrating  by  parts  the  integral  in  (A-  12).  A.ny  assumption  which 

gj 

then  expresses  7-7  in  terms  of  J  is  merely  equivalent  to  an  interpolation 
ox 

formula  for  evaluating  this  integral . 

For  these  reasons  v\e  have  found  convenient  to  use  the  direct  approach 
and  the  only  task  which  now  remains  is  tt  extern!  it  to  a  non-planar  surface  ir 
sheet.  In  place  of  (A.  7)  the  integral  ci  concern  to  us  is 


7  ~  \  J(s,)(n'  •.^H^CkrteKk::') 
k  *  1 


(A..  18) 


and  a  typical  integral  equation  is 
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(A.  19) 


'V>  =  YoK(s,J(s)  +  1  £“5c  1 

(see  eq.  4.  15). 

If  £  is  not  on  C.  we  can  interchange  the  derivative  and  integral  operations 
in  (A.  18)  to  obtain 


I  -  f 


i^jH^Qcr)  ~  (n’  •  r)  +  (n’  •  r)  H^(kr)l  d(kc’) 


JC 


Ha)(kr) 


C 

\  JtsO  •{  (#*-rKs  -9)  — “ -  +  (n'  •  r)(n  • .  ?)h'  ;  (kr)id(ks') 

Jc  L 


which  becomes 


r 

\ 


\ 

Jc 


I  =  \  J(s')(s;  •  r)(s  •r)Hil'(krNi^',) 


0 


J(s')  \  (n*  •  r)6t  •  r)  -  ts'  •  r)(s  •  r)j  (kr)d(ks') 


J, 


(A.  20) 


M)s 


on  using  the  differential  equation  for  (kr). 

To  find  the  limit  as  the  observation  point  £  approaches  C,  it  is  necessary 
to  treat  analytically  the  contribution  of  the  self  cell  fo  the  second  integral  in 
eq.  (A. 20).  Since 


lim 


.  J(s')  l(n'  •  r)(n  •  r)  -  (s'  •  r)(s  •  r)l  H  (krfd(ks')  =  +  ——  J(s) 

£-»C  \  L  11  7Tk6 


A 


where  26  is  the  width  of  the  self  cell  centered  on  s'  =  s,  we  have 
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Equation  (A.  22)  is  the  analogue  of  (A.  14)  and  (A.  16)  combined,  and  is 
the  desired  result.  Though  the  integrand  of  the  second  integral  is  a  rapidly  vary¬ 
ing  function  of  s'  for  small  s'-s,  the  particular  advantage  of  this  form  is  that 
if  J(s')  is  constant  over  a  cell  (as  is  assumed  in  the  moment  method),  the  contri¬ 
bution  of  that  cell  is  given  precisely  by  the  difference  in  the  end  values  of 
(s  •  r)H^(kr).  For  the  first  integral  in  (A.  22),  the  self  cell  contribution  vanishes 
as  the  cell  size  tends  to  zero.  The  analogy  with  the  results  for  a  planar  sheet  is 
now  complete.. 
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APPENDIX  B 


PROGRAM  RAMD 


Program  RAMD  solves  the  integral  equations  2.  9  and  2. 16  for  a  cylindrical 
body  satisfying  a  surface  impedance  boundary  condition.  Although  its  predecessor 
progiam,  RAM  IB,  accomplishes  this  same  task,  RAMD  is  a  significantly  more 
compact  and  more  efficient  program.  This  compression  was  ac;hie\cd  by  a 
reduction  in  the  number  of  arrays  created  by  the  program,  by  the  exploitation  of 
the  duality  between  the  E-  and  H-polarization  integral  equations,  and  by  the 
elimination  of  options  that  were  seldom  exercised  in  the  former  program.  As  a 
result,  RAMD  occupies  less  than  1/  3  the  core  required  by  RAM  IB,  implying  that 
the  number  of  surface  sampling  points  could  probably  be  increased  by  50  percent 
or  more.  The  program  is  manifestly  shorter  because  it  solves  the  simpler  of  the 
two  available  equations  and,  not  least,  it  produces  the  correct  solution  for  E- 
polarization  for  edged  structures  'hereas  its  predecessor  did  not. 

In  addition  to  the  MAIN  program,  RAMD  has  five  subroutines:  GEOM, 

HANK,  ADAM,  FLIP  and  ZFUN.  The  MAIN  program  reads  control  information 
from  the  input  data  stream,  creates  the  matrix  elements  for  the  geometry  at  hand, 
sums  the  surface  currents  to  obtair  the  far  scattered  fields,  and  indexes  through 
the  desired  angles  of  incidence  and  scattering.  MAIN  calls  subroutine  GEOM  to 
create  the  coordinates  of  sampling  points  on  the  body  profile  if  the  profile  is 
describable  in  terms  of  straight  or  circular  arc  segments  but,  if  necessary,  the 
sampling  points  may  be  fed  in  one  at  a  time  by  by-passing  GEOM.  GEOM  assigns 
a  specific  impedance  to  each  surface  sampling  point  according  to  two  options  and, 
should  neither  be  sil  able  to  the  user’s  needs,  a  third  option  invokes  a  call  to  sub¬ 
routine  ZFUN.  In  the  program  listing  below,  ZFUN  is  merely  a  dummy  subroutine 
necessary  for  the  successful  compilation  of  RAMD  and  should  be  replaced  by  the 
appropriate  version  required  by  the  user. 

Subroutine  HANK  generates  Hankel  functions  of  order  0  and  1  as  required 
by  the  integral  equations.  This  particular  subroutine  differs  from  previous  ver¬ 
sions  (in  REST  and  RAM  IB)  in  that  polynomial  approximations  ef  the  Hankel 
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functions  are  employed,  a  distinct  advantage  for  large  arguments;  subroutine 
ADAM  assists  HANK  by  summing  the  terms  in  the  polynomial.  Once  the  matrix 
elements  have  been  filled  in  by  the  MAIN  program,  subroutine  FLIP  is  called  to 
invert  the  matrix.  FLIP,  like  its  cousin  ZV08  in  previous  programs,  is  essentially 
a  copy  of  subroutine  MINY  from  IBM’s  Scientific  Subroutine  Package.  FLIP  also 
computes  the  surface  currents  on  the  profile  using  the  inverted  matrix.  MAIN 
then  sums  the  currents  and,  depending  on  the  particular  options  exercised  on  input, 
produces  either  a  bistatic  or  backscattering  pattern  of  the  obstacle,  for  either  E- 
or  H- polarization.. 


#*«#**#«  *  $  *  *  *  *  £  $  S  £  $  #$$$#£ 

r  INPUT  FORMAT  FOR  PROGRAM  R AMO - REVISION  OF  MAY,  197B  C 

C$*fc*$*fc«**$*#S^S*$S*S*S*  **«*£*$ 

C  CARO  1  FORMAT  I1BA4)  TITLF  CARO;  USE  IIP  TO  1?  COLUMNS  C 

C  CARO  7  FORMAT  I  I  2,  IB , F 1 0. ft  )  MORE ,  KOOF ,  7  F  AC.  C 

C  MOR  F  =0  THIS  WILL  RF  THE  LAST  RUN  FOR  THIS  OAT  A  S  FT  0 

C  MOR  F= 1  THERE  ARE  MORE  HAT  A  TO  RF  R  F  AO  AFTER  THIS  SFT  C 

C  KO0E=0  COMPHTFS  R1SIATIC  SCATTERING  PATTFRN  C 

C  K00F=1  COMPHTFS  HACKSCATTFR ING  PATTFRN  C 

C.  Z  F  AC.  A  FACTOR  MULTIPLYING  ALL  F  |  FMFNT  IMPEDANCES  C 

i  -f  $*«$*$**  I***  $  $  $  sj!  $$  *  $  $  »  *  Xc  X!  If  if  if  #  C 

C  CARO  B  FORMAT  U?,IB,E10.5)  LL,M, WAVE  C. 

C  LL  TOTAL  NHMRRR  OF  S F G M F M T S  ON  THF  PROFILE  C 

C  M  TnTAL  NUMPFR  OF  POINTS  ON  THE  PR  OF  J  L  F  C 

C  WAVE  WAVELENGTH  C. 

C  CARO  4  FORMAT  (  I ? ,  IB , ’ F 1 0. ft )  I  P P , I  OP T , F  I R S T , I  AS T , I NK  C 


IPP=1 
!  PP=  ? 

I  OPT  =0 
I OPT= 1 
FIRST 
LAST 
INK 


F-POLAR I  ?  AT  ION 
H-POL  AR I  7  AT  I  ON 

BODY  GEOMETRY  READ  IN  POINT  R  Y  POINT 
BODY  GEOMETRY  GENERATED  INTERNALLY 
INITIAL  SCATTERING  AND  INCIDENCF  AMGlF 
FINAL  ANGLF 

angular  increment 


C,  CARD  5  FORMAT  (  !2»IB,ftF10.ft)  N  ,  1  MP  ,  X  A  ,  Y  A  ,  X  fi,  YR  ,  ANG  C 

C  N  NUMBER  OF  SAMPLING  POINTS  ON  THIS  SEGMENT  C 

C.  IMP  =  -1  IMPEDANCE  GIVEN  BY  US  FR -S*  IP  PL  I F  D  SUBROUTINE  C 

C.  IMP  =  0  7S(  I  )=7A+ZB*S(  I  }**?FX  C 

c  IMP=1  7S(  I  >=7A+ZR*FXPI  -7FX*S<  I)  )  C 

C.  XA,YA,XR,YR  SFGMFNT  FNDPO I  NT S  C 

C  ANG  ANGLF  SHBTFNDtD  BY  THF  SFGMFNT  C 

C**<:*#**«*##^##S5*#*#*#*^«*«********#^*****##4:i*#**^**<t**St##*^*n:*<:#**#^*C 
C.  CARD  ft  FORMAT  IftX,ftF}0.ft)  7A,ZR.ZFX  C 

C  7  A ,  7  B  COMPLEX  IMPEDANCE  CONSTANTS  C. 

C  ZEX  REAL  IMP F n ANC F  CONSTANT  C 

r:»«*«$**^##**«^***,Jt  **#*#***#*«#«##<: #*#*>?<!*« 

C,  CARD  7  FORMAT  (I?)  INTEGER  ZERO  IN  COLUMN  ?:  SHUTS  C 

C  OFF  READING  OF  SEGMENT  PARAMETERS  C 

CysSXs***^#****#’?*  «£*$£**:)<>*$**«  X=  $«&!!:*$$>*=!':>(:  4  Xt$*S$$C. 

C  CARD  8  FORMAT  <I?»IB,F]0.ft)  MORF ,  KODF- ,  7F  AC.  C 

C  THIS  CARD  is  USFD  ONLY  IF,  UN  CARD  ?,  M0RF=1  C 

CXtSWJfeXt*#*^#^*#***#*#**#***#**^#******^*****#:*#:!:*!*"  *$£*«*$*$***  *!><si!;:!;?*!#$$C 
C  CARD  9  FORMAT  (  I?,IB,4Fl0.ft)  LHM  p  (  I  ,  1)  ,  L  UM  °  (  I  ,  7  )  .  X  (  I  )  ,  Y  (  I)  ,  7  S  (  I)  C 
C.  LUMP  (1,1)  CFLL  ID  NUMBER  C 

C  LUMP (I,?)  SFGMFNT  ID  NUMBER  C 

C  XII>,Y(I)  CELL  COORDINATES  C 

C  ZS(I)  COMPLFX  IMPFDANCF  of  THF  CFLL  C 

C  NOTE:  CARD  9  IS  USED  ONLY  JE  IOPT=0  AND  THFRF  MUST  RF  DNF  SUCH  CARD  C 
C  FOR  FVFRY  CFLL  ON  THE  BODY.  SIMILARLY,  CARDS  ft, ft  AMD  7  ARE  USED  C 
C  ONLY  IF  I OPT= 1  AND  THFRF  MUST  BE  ONE  FACH  OF  CARDS  ft  AND  ft  FOR  EACH  C 
C  SEGMENT  SPECIFIED  ON  INPUT,  C 
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r omp i  f x  a n on, ion ,phi  ( ion)  ,pink(  i nos^uM, dfl ,ri  ,R2» 7S(  ino) 
RFAL  LAST, INK 

niMFNS  ION  X  (  100  >  ,  Y(  ’.no  )  ,XN(  1  00)  ,  YN(  100)  ,S  UOO  )  ,DSQ(  100  ) 

0 1  MFNS I  DM  I  Od  R),)  I  IMP  <  ion,?),  IP(U  (?) 

DATA  RFn.nin, I  POL /O.Ol 745 T?Q, 57.2957R, • Ff FF  •  ,  1 HHHH 1 / 

C . R  f  AO  INPUT  OAT  A  AMO  GFNFRATF  BODY  PROF  I L  F 

s  r f ao  <5,1  no:  io 

R  F AO  (5,200)  MOPF, KOOF, 7FAC 

RFAO  (5,2 00  )  LI  ,M,WAVF 

R  F  AO  (5,200)  IPP, IOPT, FIRST, LAST, INK 

XK=6.?RS! R5/WAVF 

WR I T  r  (6,150)  10 

I F  ( IPPT.FO.O )  GO  TO  10 

WR i T  F  (6, TOO) 

CALL.  OFDM  U 'IMP,  X,Y,XN.,YN,<  .  OSO ,  ?  S  ,  M  .  L  L  ) 

GO  TO  20 

10  RFAO  (5,200)  LUMP ( 1 , 1 ),l"MP( 1,2 ) tX( 1 ) ,Y(1 ),ZS( l ) 

1  )  =  ;0  .  0 
on  15  1=2, M 

RFAO  (5.200)  (L'>MP(  I,  J),  J=1,2).X(  I),Y(  I),ZS(  I  ) 

t x = x (  n-x ( i-i ) 

TY-YI I )-Y (1-1 ) 

0S0(  I  )  sSORT  <  TX*TX«-TY*  l\  ) 

s ( i )  =  s ( i - 1 )  +  nso<  i ) 

XW(  I  )=-TY/OSO(  I  l 
16  YN( I >=TX/OSO( I ) 

20  rF  (KOOF.NF.O)  00  TO  25 
NINOl 

NRIT  =  1  +  IFIX((I.AST-F1RST)/INK) 

on  to  to 

2  5  NR  I T  =  0 

NINOl  +  l  FI  Y(  (I.  AST-F  IRST  I  /  INK  1 
TO  WRITF  (6,«f.0)  IPOL(  IPP),LL,M,NINC,NRIT,WAVF 
on  T5  I=1,M 
55  OSO(  I  )=0S0(  I  ) / W  AV  F 

r, . CONSTRUCT  MATRIX  F|  FMFNTS 

TB  00  55  1=1  ,M 
nn  55  .tsi, m 
IF  ( ! . rO, J )  on  TO  AO 
T  X  =  X  (  I)-X(  ) 

TY  =  Y (  I  )  —  Y ( J  ) 

P * SDK  T ( TX#TX4TY*  1 Y  ) 

RPO=P*XK 

CNR  =  -(  TX*XM(  ,1  l+TY*  YM  .))  )  /P 
CALL  HANK  (RPQ,  1  ,R.I,RV  ) 

R 1  =  1 . 5 70 706* PS 0(  J)*C.Nk^C,MP|  X(-BY,RJ) 

CALL  HANK (RPO.O.R J,RY ) 

R2=1.570796*0S0(  J  )  *('. MPL  X  (  R.l ,  R  Y  ) 
on  TO  46 


40  R1=CMPLX<0.5,0.0) 

B2=r>SQ<  J  )*CMPLX(  1.  5 70794, 0.0? 8 7985+ A L0G<  DS0(  J  )  )  ) 

44  I  F  (  1PP.FQ.1  )  on  TO  BO 
AU,J)  =  R1+R?*7S(J)*7FAC 
on  to  55 

50  A( I , J)=R?+R1*7S( J)*7FAC 
55  CONTINUE 

C . COMPUTF  INCIOFNT  F  I  FLO  AMO  INVERT  MATRIX 

tfta=rfo*first 
CT=COS( TFT  A  ) 

5 T  =  S I M < TPTA  ) 

nn  40  i = i , m 

HOLOs-XK* (CT*X ( I)+ST*Y(in 
40  P I NK (  I  ) =CMPL X( cns< HOLD), SI N< HOLD)  ) 

CALL  FLIF(A,M. PINK, PHI, II 

C . PRINT  OUT  CURRENTS  AND  FLFMFNT  PkOPFRTIFS  FOR  FIRST  ANGLE 

WRITF  (4,150)  10 

WRITF  (4,400) 
nn  45  1  =  1, M 
AMP=CARS( PHI ( I ) ) 

PHASF  =  n  I  G*  ATAN?  (  A  I M  AC,  (  Ph  I  (  f  )  )  , R F  AL  (  PHI  (  I  )  )  ) 
nFL  =  ZFAC*Z.S(  I  ) 

45  WRITF  (4,750)  (  LUMP  (  I  ,  J  )  ,  J=  1 ,  ?  )  ,  X  (  I  )  ,  Y  (  I  )  ,  S  (  I  ) ,  OSO  (  I  )  ,  AMP ,  PHA.SF  , 

F.nFL 

C . nnPF  nnT  thf  appropriate  fiflo  factors 

THF=FIRST- I MK 
IF  (KOOF.FG.l )  GO  TO  70 
WRITF  (4,800)  F I  PS  T 
GO  T 0  75 

70  WRITF  (4,400) 

75  THF  =  THF-*- 1  NK 

IF  ( THF. GT. LAST )  GO  Tn  105 
IF  ( THF,  FO.  FIRST  )  GO  TO  RF 
TFTA=RFn*THF 
CT=CnS ( TFTA ) 

ST  =  S I N ( TFTA) 

C . IN  THF  FOLLOWING  LOOP,  PINK  IS  NOT  NECFSSARILY  THF  INCIOFNT  FIFLO 

no  RO  J=l,M 

HOLn=-XKS  (  C,T*  X  (  J)+ST*Y(.))  ) 

80  P  I  NK  (  J)=C.MPLX(CnS(  HOL 0  )  ,  S  I  N  (  HnL  0  )  ) 

IF  (KOOF.FO.O)  GO  TO  R5 
CALL  FLIP( A, M, PINK, PHI,?) 

85  SUM=CMPLX(0. 0,0.0) 

c . Ann  op  thf  currents 

nn  95  j=i ,m 

of  l  =  nso  (  ,1  )  *  PHI  (  J  )*  PINK  (  vl  ) 

IF  ( I PP .FO. 1 )  GO  TO  90 

S(IM=SUM+OFL*  (ZFAC*7S(.I ) -C T*XN  (  J  ) -ST *YN  ( ,1  )  ) 

GO  TO  94 

90  SUM=  SUM  +  OF  L# ( 1 . 0-7 F AC* 7S ( J ) * ( CT~ XN( J ) +ST* YN( J ) ) ) 

95  CONTINUE 
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scat  =  ?o.o*ai.ogio(cars  (  siimii  +  1 .9612 

WRITF  (<S,900)  THF.  SCAT 

an  to 

ins  if  imorf.fq.o)  an  to  s 

RFAD  (S,?00)  MORF ,KOnF,7FAf 

an  m  tb 

inn  FORMAT  (1RA4) 

150  FORMAT  (lHltlSAM 

?oo  format  ( i?,  n,SFin.s> 

?so  format  (?i s, sFin.s,Fin.R,?Fin.s) 

500  FORMAT  (10HOSFG  NilM,llX,?4HFNI)P0(NT$  OF  THF  SFGMFNT  ,  1 1 X  , 

F,  1 RHS  FGMFMT  PARAMFT  FRS/ 1 1 H  MUM  C  F  L  L  S  ,  <5X  ,  ?HX  A  ,  ft  X  ,  PHY  A  ,  8  X  ,  ?HX  R  t  8  X  , 
£?HYR»ftX,21HAMGl  F  RAO II IS  I  FNGTH/ ) 

400  FORMAT  ( //?! X- 14HKFY  PARAMFTFRS// 

F.  1  AX  *?  1  H  I  Nf.  I  DFNT  PHI.  AR  I  Z  AT  1  ON  ,  22  X  ,  1  A  1  /  / 
f.  1 6  X  »  2  RH  N*  IM  R  F  R  OF  SFGMFN1S  IISF0«I?1// 

F.1F>X,^4HT0TAL  NIIMRFR  OF  POINTS  ON  THF  ROPY,  I  10// 

F.16X,  35HNIIMRFR  OF  INC.IOFNT  FIFLD  0 1 R FCT I ONS , l 9/ / 

F.  1  6X  ,  ?9KNI  IMR  FR  OF  PISTATIC  0 1 R FCT I ONS » I  1 5/ / 

F.1F.X,  10HWAVFI  FNGT h«  FB4.S) 

SOP  FORMAT  {  1 1  HO  I  SF  G,  4X  ,  4HX  (  I)  ,  F>X  ,  4HY  (I)  t  F>X  ,  4HS  I  II  ,  5  X  ,  6HDS0  <  II  , 
F.4X  ,  SHMOOI  J  ),  4X,  SHARP,!  J  )  ,4X,  SHRSI  I)  ♦  5  X  ,  SH  XS  (I  )  /  ) 

(Son  FORMAT  (  1H1  ,??X,?RHRACKSC  ATTFRINO  CROSS  SFC  T  I  0N//24X  , 

F.pfiHTHFTA  10*LOG<  S  I  GM  A/L  AMRO  A  )  /  ) 

ROO  FORMAT  (  1 H 1  , 1  9  X  ,  5  SHR  I  ST  AT  I C,  SCATTFRING  CROSS  SFCTI0N/1PX, 

F.29HFOR  INCIOFNT  FIFlO  OJRFCT  I0N=»F7,2//?4X, 

F.pHHTHFTA  lO^LOOI  SIGMA /L  AMRO  A)  /  ) 

900  FORMAT  <  lftX,Fn.2,FlS,?) 

FMO 


SURROUTINF  GFOM( LUMP,X,Y,XN,YN, S»OSO,ZS,M,LL  ) 

CHMPLFX  Z  S  ( 100 ) . Z  A, ?R 

01 MFNSION  X(  100)  ,  Y(  100  )  ,  XN(  100  )  ,  YN(  1  00  )  ,  0S0<  100)  ,  S<  100) 

OIMFNS ION  LUMP! 100,2  ) 

OATA  RFn/0.0174S32P/ 

1=0 
L  =  0 

C . R  F AD  INPUT  PAR AMF  TFRS  ANO  PRFPARF  TO  GFNFRATF  SAMPLING  POINTS 

10  R  F  AD  (5, POO)  N, IMP, XA, YA,XB, YR, ANG 
IF  (N.FO.O)  GO  TO  IPO 
LIM=P*N-1 

R FAR  (  S , 2  SO  )  ZA,ZR,ZFX 
TX  =X  R-X  A 
TY=YR— YA 

0=SQRT(TX*TX+TY*TY > 

IF  (ANG. FO, 0.0)  GO  TO  ?0 
T  =  0. 5*R  FO*  ANG 
TRX=TX+TY*COTAN(T  ) 

TRY=TY-TX*COTAN( T ) 

RA0=0.S*0/S  IN( T  ) 

ARC  =  2. 0*R AO*  T 

alf=t/n 

oi o=2.o*rao* alf 

GO  TO  30 
20  R  A0=999. 999 
ARC  =  0 

ni o=o/n 

C . START  GFNFRATING 

30  00  UO  J  I  M=  l ,  2 

L  =  L  +  1 

no  ioo  j=i,lim,p 

1  =  1  +  1 

LUMP) 1 , 1  )  =  I 
LUMP ( I ,2  )=L 

IF  (I.F0.100)  WRITF  (6,400) 

IF  (JIM.F0.2)  GO  TO  PO 
IF  <  ANG. FO. 0.0)  GO  TO  40 
SINO=SIN( J*ALF  ) 

cnso=cns{ j*alf  ) 

X ( I )=XA+O.S*<  TRX* ( 1 .0-C0S0)-TRY*SIN0) 

Y ( I )=YA+O.S*(TRX*SING+TRY*( l.O-COSQ) ) 

XN( I ) =-0 . S* (TRX*COSO+TRY*SIN0)/RAD 
YN ( I  )  =  O.S*{ TRX*SINQ-TRY*CnSO)/RAO 
GO  TO  SO 

40  X(  I  )  =  XA  +  0. S* J*TX/N 
Y ( I  )=YA+O.S*J*TY/N 
XN<  I  )=-TY/U 
YN( I ) =  TX/O 

so  s  ( i ) =0.  s*  j*o  i  n 
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C . CDMPIITF  THF  IMPFOANCFS 

IF  (  IMP  )  (SO, 70,  80 

60  CALL  ZFUN( ?A,ZR, 7FX, S(  I)  ,ZS(I)  } 
on  m  ion 

70  z  S  (  I  )  =7  A  +  Z  B*  S  (  I  }  *  *  Z  F  X 
on  to  ioo 

80  7S<  I  >=ZA+ZB*FXP(-7FX*S<  5  )  ) 

on  m  ioo 

C . FROM  H  CR  F  TO  100  WF  CRFATF  THF  SEGMENT  IMAOF 

90  K  =  I  -N 

xm=x(K» 

Y(  I  )=-Y(K  ) 

X  N  (  I  )  =  X  N  (  K  ) 

YN(  i )=- YNIK  ) 

s  ( i  >  =  s ;  k  ) 

ZSI  I  ) =7 F ( K  ) 
ioo  nsQ( I i - n i n 

if  (  jim.fo.i  )  on  in  no 

Y  A  =  -Y  A 
YR=-YR 

110  WRITF  (6,300)  L,N, XA, YA,XR, YB, ANG,RAD,ARC 
on  TO  10 
1 pn  m= i 
LL  =  L 

POO  FORMAT  (I?, 13.SF10.S) 

PRO  FORMAT  (SX,RF10.R) 

300  FORMAT  (  I  3,  16, 3X,4F10- S,FR.P,F8. 3,F8.4  ) 

600  FORMAT  ( 36H0WARN I  NO :  8F'VF  GFNFKATFO  100  POINTS/) 
RETURN 
FNO 
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SURROilT  1  NF  HANK  (R  ,  N  ♦  R  .1 »  R  Y  ) 

C . SHRROtlT  1  NF  REQUIRES  R>0  ANn  N  EITHER  0  HR  1 

DIMENSION  A(M*R(7)*CI7)tD*7)*E(7l«F(7)»(S(7)«H(7) 

DATA  A t P, C , D, Ft F, 0 , H/ 1 .0,-?. 2499997, 1.2656208, -0.31 63866, 
£0.0444479,-0.0039444,0.00021 ,0.36746691 ,0 . 60559 366 , -0 . 74350384 , 
£0.  253 00 117, -0.04 261 214,0.004 27916, -0.000 24 84 6,  0.  5,-0.  8249985, 
£0.21093573,-0.03954289,0.00443319,-0.00031781,-0.00001109, 
£-0.8368198,0.2212091,2. 1882709,-1.3184827,0.3123951,-0.0400976, 
£0.0027R7 3,0.797 88456, -0.00000077 ,-0.0055274, -0.0000951 2, 

£0. 00137237,-0. 00072805,0.00014476,-0.78539816,-0.04166397, 

£-0. 000039 R4, 0.002 62 57 3, -0.00054 12 8, -0.00029 333, 0.0001 3558, 
£0.79788486,0.00000156,0.01659667,0.00017105,-0.00249511, 

£0.00 11 3653, -0.000 20033 ,-2.356 19449,0. 12499612,0.0000565, 
£-0.00637879,0. 00074348,0.00079824,-0.00029166/ 

IF  (R. IF. 0.0)  00  TO  50 
IF  (R. OT. 3.0)  on  TO  20 
X=R*R/9.0 

IF  (N.MF.O)  on  TO  10 
CALL  ADAM ( A ,X ,R J ) 

CALL  ADAM ( R, X , Y 1 

RY  =  0.636619 8*  AL no  (0,5*R  ) *R.J+Y 

RETURN 

10  IF  <N.NF.l )  OH  TO  50 
CALL  ADAM ( C, X, Y ) 

R J=R*Y 

CALL  ADAM ( H, X, Y ) 

RY=0. 636619 R^AL DO (0.5*R )*BJ+Y/R 
RETURN 
20  X=3.0/R 

IF  (N.NF.O)  00  TO  30 
CALL  ADAM  (  F  , >'  ,  Y  ) 

FOOL  =Y /SOR T ( R  ) 

CALL  ADAM ( F . X , Y  ) 

On  TO  40 

30  IF  (N.NF.l )  On  Tn  50 
CALL  ADAM( 0, X,Y) 

FnnL=Y /SORT ( R ) 

CALI  ADAM( H. X, Y ) 

40  T=k+Y 

Rj  =  Fnni *cns( T ) 

RYsFOOL- S  I N  (  T  ) 

RFTIIRN 
50  N= 1 00 
RFTIIRN 
END 
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f 

f 

1 

- 

I 

1 

r- 

p 

jfe 

it 

M 

StlRROIITlMF  FLIP!  A,N,X,Y,  1  AT  ) 

i k 

i 

CHMPLFX  A( 100, 101 ),X( 100)»Y( 1 00 ) t 0 , R I 0 A, HOLD 

6 

0 1 MFNS I ON  L ( 100  1  ,M( 100 ) 

1 

IF  ( I AT.GT.l )  00  TO  130 

p 

p=C  MPL  X I 1 .0,0.0 ) 
nn  ho  k=i.n 

fc 

L«K)=K 

1 

M( K  )  =  K 

B1GA=A(K.K) 

i 

nn  ?o  j=k,n 

W/  } 

nn  ?o  i =k  . n 

10 

IF  (CARS(Rir,A).OF.CARS(  A(  I  ,.l )  1  )  0,n  fO  ?0 

[ 

fl 1 0 A  =  A ( I  ,  J) 

tk 

L  (  K  )  =  1 

i 

M(K  >  =  J 

1 

?o 

CnNTINIlF 

J  =  L  IK  ) 

* 

if  ij.lf.k)  r,n  Tn  33 
nn  30  1  =  1 ,  n 

Hni0=-A ( K,  I  ) 

i 

AIK,  I  )  =  A(.)t  I  ) 

I 

30 

A(  .1,  I  )=HOL  n 

I 

33 

1 =  M(K  ) 

IF  I  I  ,  LF.K  )  0,n  TO  A3 

I 

nn  40  j=i  ,n 

HOI n  =  -A ( J, K  ) 

1 

A  I  J  ,  K  )=A(.i,  I  ) 

I 

40 

A( J, I ) =HOLn 

■ 

4  3 

IF  (CARS  I  RIO.A  )  ,NF  .0.0  )  0,0  TO  30 

1 

r:=CMPixio.  o.o.o) 

1 

* 

i 

RFTURN 

1 

i 

30 

nn  33  1=1, n 

P 

1 

| 

IF  I  I  .FQ.K  )  0,n  Tn  33 

S 

A  (  I  ,  K  )  =- A  I  I  ,K  >  /RIO,A 

§ 

33 

OONT INUF 

K 

nn  33  i=i,n 

is 

nn  33  „i=i, n 

I 

IF  (  I.FO.K.nR.  J.FO.K  )  0,0  TO  33 

it 

AI  I  ,  J  )  =  A  I  I,K)*A(K«.!>  +  A(  I  ,  J  ) 

i 

6  S 

CONTINUF 

L# 

{ 

00  73  .1  =  1, N 

g* 

IF  1  J.FO.K  )  0,0  TO  73 

|  IP 

%  £r 

73 

A(K,J)=A!K,J)/RIGA 

CONTINUF 

1  ’* 

! 

n=n*Rir,A 

If 

i 

,? 

HO 

A(K,K  )  =  1.0 /RIO, A 

II 

t  -  > 

Eh 

1  V, 

1  v. 

100 

K  =  K-1 

K  i 

i  £ 

IF  (K.LF.O)  0,0  TO  130 

all 

’  -S'.. 

I  =  L  I  K  1 

1 

i  5r- 

i  S* 

IF  II  .LF.K  I  0,0  TO  1  ?0 

1 

!  S 

03 

i 

nn  no  j  =  i,n 
HOLn= A ( J,  K  ) 

A  <  J  »  K ) =-A ( J  » !  ) 
lio  A(j,n=Hnin 
120  J  =  M ( K  ) 

if  <  J. lf.k  )  r,n  rn  ion 
on  130  1  =  1, N 
HOL 0=  A ( K , I  ) 

A ( K , I  )  =-A ( J  ,  I  ) 

1 30  A(j,n=Hni_n 
on  m  ion 
iso  nn  200  i=i, n 

Y ( I )=CMPLX(0. 0,0.0) 
on  200  J=  1 ,  N 

200  V(I)*AII,J|*X(J)  +  Ym 
RFTURN 
FNO 
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SUBRni  IT  I  NF  ADAM  ( C  *  X  «  Y  ) 

ni MFNSI HN  C( 7 ) 

Y  =  X*C.(  7  ) 

nn  in  i=i,5 

in  y=x* ( r. ( 7  - 1 )  +y  ) 

y=Y+r.  ( i ) 

RFTURN 

FND 


SIIRROUT  I  MF  7  Fl IN'  7A,  ZR»  ZFX,  7.S  ) 

7  WA+ZR  +  ZFX 

RFTIIRN 

FND 


REFERENCES 


Andreasen,  M.G.:  (1965)  ’’Scattering  from  cylinders  with  arbitrary  surface 
impedance"..  Proc.  IEEE’  53,  812-817. 

Bowman,  J.J.  (1967),  "High-frequency  backscattering  from  an  absorbing  infinite 
strip  with  arbitrary  face  impedarces",  Can,  J,  Phys.  45,  2409-2430, 

Knott,  E.F.,  V.V.  Lie  pa  and  T,.  B.  A,  Senior  (1973),  "Non-specular  Radar  Cross 
Section  Study",  The  University  of  Michigan  Radiation  Laboratory  Report 
No.  011062-1-F  (AFAL-TR-73-70),  April  1973. 

Knott,  E.  F,  and  T.  B.  A.  Senior  (1973),  "Non-specular  Radar  Cross  Section  Study", 
The  University  of  Michigan  Radiation  Laboratory  Report  No.  011062-1-T 
(AFAL-TR-73-2),  February  1973. 

Laxpati,  S  R,  (1973),  "Derivation  of  integral  equations  for  a  two  dimensional 
scatterer  with  impedance  boundary  condition  and  partially  clad  by  an 
absorptive  sheet".  The  University  of  Michigan  Radiation  Laboratory 
Memorandum  No.  0117  64-  505-M. 

Lebedev,  N.  N.  and  I.  P.  Skal'skaya  (1963),  "A  new  method  for  solving  problems 

of  diffraction  of  electromagnetic  waves  by  a  wedge  with  finite  conductivity", 
Sov.  Phys.- Tech.  Phys.  867-874 

Maliuzhinets,  G.D.  (1960),  "Das  Sommerfeldsche  Integral  und  die  Losung  von 
Beugungsaufgaben  in  Winkelgebieten",  Ann,  Physiki),  107-112, 

Oshiro,  F.K.,  K.M.  Mitzner,  S.  S.  Locus,  et  al,  ( 1 97 1) .  "Calculations  of  Radar 
Cross-section",  Northrop  Corporation,  Norair  Division  Technical  Report 
No.  AFAL-TR-7  1-50. 

Senior,  T.B.A.  (1952),  "Diffraction  by  a  semi-infinite  metallic  sheet",  Proc.  Roy.. 
Soc.  213A,  436-458. 

Senior,  T.B.A.  (1959),  "Diffraction  bv  an  imperfectly  conducting  wedge",  Commun, 
Pure  Appl.  Math,  J 2,  337-372, 

Senior,  T.  B.  A.  (1960),  "Diffraction  by  an  imperfectly  conducting  half-plane  at 
oblique  incidence",  Appl.  Scl.  Res.  ^B,  35-61. 

Senior,  T.B.A.  (1962),  "A  note  on  impedance  boundary  conditions".  Can.  J. 

Phys.  40,  663-665. 


102 


REFERENCES 


Andreasen,  M.G.  (1965)  ''Scattering  from  cylinders  with  arbitrary  surface 
impedance".  Proc.  IEEE  _33,  812-817. 

Bowman,  J.J.  (1967),  "High-frequency  backscattering  from  an  absorbing  infinite 
strip  with  arbitrary  face  impedarces",  Can.  J.  Phys.  45,  2409-2430. 

Knott,  ELF.,  V.V.  Liepa  and  T.  B.  A.  Senior  (1973),  "Non-specular  Radar  Cross 
Section  Study",  The  University  of  Michigan  Radiation  Laboratory  Report 
No.  011062-1-F  (A FA L-T R-7 3-70) ,  April  1973. 

Knott,  E.  F.;  and  T.  B.  A.  Senior  (1973),  "Non-specular  Radar  Cross  Section  Study", 
The  University  of  Michigan  Radiation  Laboratory  Report  No.  011062-1-T 
(AFAL-TR-73-2),  February  1973. 

Laxpati,  S  R..  (1973),.  "Derivation  of  integral  equations  for  a  two  dimensional 
scatterer  with  impedance  boundary  condition  and  partially  clad  by  an 
absorptive  sheet".  The  University  of  Michigan  Radiation  Laboratory 
Memorandum  No.  0117  64-  505-M. 

Lebedev,  N.  N.  and  I.  P.  Skal'skaya  (1963),  "A  new  method  for  solving  problems 

of  diffraction  of  electromagnetic  waves  by  a  wedge  with  finite  conductivity", 
Sov.  Phys..  -  Tech.  Phys.  7,  867-874 

Maliuzhinets,  G.  D.  (1960),  "Das  Sommerfeldsche  Integral  und  die  Losung  von 
Beugungsaufgaben  in  Winkelgebieten",  Ann.:  Physik  6,  107-112. 

Oshiro,  F.K.,  K.M.,  Mitzner,  S.  S.  Locus,  et  al.  (i971),  "Calculations  of  Radar 
Cross-section",  Northrop  Corporation,  Norair  Division  Technical  Report 
No.  AFAL-TR-7  1-50. 

Senior,  T.B.A.  (1952),  "Diffraction  by  a  semi-infinite  metallic  sheet",  Proc.  Roy. 
Soc,  213A,  436-458. 

Senior,  T.B.A.  (1959),  "Diffraction  by  an  imperfectly  conducting  wedge",  Commun.; 
Pure  Appl.  Math.  _L2,  337-372. 

Senior,  T.  B.A.  (1960),  "Diffraction  by  an  imperfectly  conducting  half-plane  at 
oblique  incidence",  Appl.  Scl.  Res.  oB,  35-61. 

Senior,  T.B.A.  (1962),  "A  note  on  impedance  boundary  conditions".  Can.  J. 

Phys.  40,  663-665. 


102 


UNCLASSIFIED 


fteiuM’t  Ci«i*Mlire*tion 


Security  Claamification 


